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The hallmark of superconductivity is the diamagnetic response to external magnetic field. In strik- 
ing contrast to this behavior, a paramagnetic response or paramagnetic Meissner effect was observed 
in ceramic high-Tc and in conventional superconductors. The present review is given on this inter- 
esting effect and related phenomena. We begin with a detailed discussion of experimental results 
on the paramagnetic Meissner effect in both granular and conventional superconductors. There are 
two main mechanisms leading to the paramagnetic response: the so called d-wave and the flux com- 
pression. In the first scenario, the Josephson critical current between two d-wave superconductors 
becomes negative or equivalently one has a tt junction. The paramagnetic signal occurs due to the 
nonzero spontaneous supercurrent circulating in a loop consisting of odd number of tt junctions. In 
addition to the d-wave mechanism we present the flux compression mechanism for the paramagnetic 
Meissner effect. The compression may be due to either an inhomogeneous superconducting tran- 
sition or flux trap inside the giant vortex state. The flux trapping which acts like a total nonzero 
spontaneous magnetic moment causes the paramagnetic signal. The anisotropic pairing scenario 
is believed to be valid for granular materials while the flux trap one can be applied to both con- 
ventional and high-Tc superconductors. The study of different phenomena by a three-dimensional 
lattice model of randomly distributed tt Josephson junctions with finite self-inductance occupies 
the main part of our review. By simulations one can show that the chiral glass phase in which 
chiralities are frozen in time and in space may occur in granular superconductors possessing d-wave 
pairing symmetry. Experimental attempts on the search for the chiral glass phase are analysed. 
Experiments on dynamical phenomena such as AC susceptibility, compensation effect, anomalous 
microwave absorption, aging effect, AC resistivity and enhancement of critical current by external 
electric fields are considered. These phenomena were studied by Monte Carlo and Langevin dy- 
namics simulations which show satisfactory agreement with experiments. We present a resistively 
shunted junction model describing rich dynamics of Josephson junction networks. 
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I. INTRODUCTION 

The phenomenon of superconductivity was first dis- 
covered in 1911 by H. Kamerhngh Onnes. He observed 

that when mercury was put into a hquid hehum bath 
and cooled down to very low temperatures, its resistiv- 
ity dropped to zero once the temperature becomes lower 
than a critical value Tc- It was soon observed that su- 
perconductors, except from zero resistivity, also exhibit 
unusual magnetic properties. In 1933 W. Meissner and 
R. Ochsenfeld found that when a superconducting speci- 
men is placed in a low enough magnetic field and is sub- 
sequently cooled through its transition temperature, the 
magnetic flux is totally expelled from its interior. This 
phenomenon was called Meissner effect. In the conven- 
tional Meissner phase the magnetic susceptibility is neg- 
ative or the response to the external magnetic field is 
diamagnetic. 

The revolutionary discovery of high temperature su- 
perconductivity by Bednorz and Muller [1] ignited an 
explosion of efforts in study of ceramic materials. One 
of such efforts led to the discovery of the so called para- 
magnetic Meissner effect (PME). Contrary to the stan- 
dard diamagnetic response to the external field, a para- 
magnetic signal was observed in certain ceramic super- 
conductors upon cooling in low enough fields [2-4]. In 
other words, the field cooled (FC) susceptibility becomes 
positive whereas the zero field cool (ZFC) magnetization 
remains negative. This effect is now referred to as the 
PME or as the Wohlleben eflPect [6]. The existence of 
the PME is quite surprising because according to the 
Ginzburg-Landau phenomenology for continuous media 
[7], the lowest free energy states should be diamagnetic. 

The anomalous PME in ceramic superconductors has 
been interpreted in the framework of different models, 
such as spontaneous supercurrents due to vortex fluc- 
tuations combined with pinning [2], orbital glass [8], the 
presence of so called tt contacts [4,9-11,6,12], and Joseph- 
son junctions [13,14]. In this review we focus on the 
most prominent d-wave mechanism proposed by Sigrist 
and Rice [9,10]. According to their theory , the nature 
of the unusual paramagnetic behavior may be related 
to the existence of tt-J unctions characterized by negative 
Josephson couplings (conventional 0-j unctions have pos- 
itive couplings). It was argued that such 7r-junctions are 
a consequence of dx2-y2 type pairing symmetry [9]. If 
a loop consists of odd number of tt junctions then one 
can show that the spontaneous current (or magnetic mo- 
ment) circulating around it becomes nonzero. Such a 
frustrated tt loop would, therefore, behave like a param- 
agnet and have the paramagnetic response to the external 
magnetic field. The existence of the d-wave pairing seems 
to be supported by phase-sensitive and phase-insensitive 
experiments [15]. 

The situation becomes ambiguous when the PME was 



observed even in conventional superconductors such as 
Nb [16,17] and Al [18]. This has prompted the appear- 
ance of a less-exotic mechanism based on flux capture 
inside a superconducting sample and its consequent com- 
pression with lowering temperature [19-21]. The flux 
trap can be caused by inhomogeneities [19,20] but could 
also be an intrinsic property of any finite-size supercon- 
ductor due to the sample boundary [21]. The trapped 
fluxes play the same role as spontaneous supercurrents 
in the d-wave mechanism and give rise to the PME. 

At present, it is not clear if one needs the exotic d- 
wave mechanism to explain the PME in ceramic materi- 
als or the flux trapping is adequate for both conventional 
low- Tc and high-Tj, superconductors [18,22]. The study 
of dynamical phenomena such as the microwave absorp- 
tion (MWA) [3], the compensation effect [23], the aging 
[24] and enhancement of critical current [25], may shed 
light on this important issue. The anomalous MWA, e.g., 
was observed [3] only in the ceramic samples showing 
the paramagnetic signal. Simulations by different groups 
show that all of dynamical effects may be captured by 
the multi-loop model of the Josephson network [12] with 
TT-junctions. So the study of dynamics effects in ceramic 
superconductors supports the rf-wave mechanism. 

Recently we have proposed that a novel thermody- 
namic phase might occur in zero external field in uncon- 
ventional superconductors [26,27]. This phase is charac- 
terized by a broken time-reversal symmetry and called a 
chiral glass phase. The key idea is that the random distri- 
bution of 0- and 7r-junctions between adjacent d^^-yz su- 
perconducting particles causes the frustration effect and 
induces the chiral glass transition. In the chiral glass 
phase, similar to spin glasses, chiralities [28] are frozen 
in space and in time. Using the multi-loop model [12] one 
can show that this phase is stable even under the screen- 
ing [27]. Since the existence of the chiral glass would 
favour the d-wave mechanism for the PME, its study is 
very important for classification of the symmetry of the 
order parameter. The experimental search for the chiral 
glass phase in ceramic samples of different groups gives 
controversial results [29-32] . In this review we also briefly 
discuss the difference between the chiral glass and other 
glassy phases in impure type-II superconductors such as 
the gauge glass, vortex and Bragg glass. 

It should be noted that results obtained before 1995 
were summarized in the review article of Sigrist and Rice 
[10]. However, these authors discussed only the rf-wave 
mechanism of the paramagnetic Meissner effect and re- 
lated experiments in ceramic materials and interesting 
dynamical phenomena correlated with the PME were not 
addressed. The present review focuses on the progress 
achieved after 1995. The main developments of this pe- 
riod are: 

a) the discovery of the PME in conventional supercon- 
ductors [18]; 

b) more detailed study of the PME and related dy- 
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namical phenomena in ceramic materials; 

c) flux compression mechanism [19,21] which is not 
based on the d-wave pairing symmetry; 

d) theoretical and experimental search for a chiral glass 
phase. 

The article is organized as follows. In section 2 we 

present a brief summary of experimental results on the 
PME in both ceramic and conventional superconductors. 
In section 3 we review the mechanisms for occurrence of tt 
jimction. One of them is relied on the scattering on mag- 
netic impurities and the another, more intrinsic, mech- 
anism is based on the anisotropic pairing of supercon- 
ducting electrons. Section 4 concerns theory and simula- 
tions of the PME in ceramic superconductors employing 
single-loop and multi-loop models of the Josephson junc- 
tion network with tt junctions. In section 5 we discuss 
the flux compression mechanism for the PME based on 
the Bean critical model and the Ginzburg-Landau equa- 
tion approach. Section 6 is focused on the simulation and 
experimental search for the chiral glass phase in granular 
materials. Experiments and simulations on dynamical 
phenomena related to the PME are presented in section 
7. The multi-loop model is shown to provide an unique 
tool to describe all of these effects. 



interior and x approaches — l/47r as T ^ 0. It is shown 
schematically in Fig. 1. Above Tc one has the paramag- 
netic behavior x = c/T [7]. The susceptibility is slightly 
positive even in the narrow temperature region below Tc 
but it has nothing to do with the PME where the FC 
susceptibility remains positive up to T = 0. 



-1/477 
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FIG. 1. Schematic temperature dependence of the suscep- 
tibility of superconductors without the PME 



II. EXPERIMENTS ON PARAMAGNETIC 
MEISSNER EFFECT 



A. High-temperature superconductors 



In the Meissner phase the total magnetic fleld B inside 
conventional superconductors vanishes and we have 



B = H + AttxM = , 
_ dM _ 1 
^ ~ 'dW ~ ' 

X arc the external field, 



(2.1) 



where H, M and x ^-re the external field, the induced 
magnetization and the susceptibility. Due to the finite 
London penetration depth (of order of 500 A) the mag- 
netic fluxes do not, however, entirely pull out from the 



Approximately 100 different cuprate materials, many 
of which are superconducting, have been discovered since 
1986. Cuprates have layered structure, with quite sim- 
ilar a- and b-ax.es and considerably longer c-axis. They 
consist of Cu02 planes, normal to c-axis, separated by 
layers containing other types of atoms. The supercon- 
ductivity occurs in the Cu02 layers and it is sensitive to 
the oxygen content. The layered structure gives rise to 
strong anisotropy of many physical quantities. 



TABLE 1. Ceramic superconductors showing the PME 









grain size 


compound 


source 


(K) (A/cm2) 




Bi2Sr2CaCu208.i82 


[4] 


84 lO'^ - 10'^ [33] 


2 - 3 


Bii.73Pbo.27Sr2Ca2Cu30y 


[4] 


110 10^ - 10^ 


4 


YBa2Cu307-5 


[34] 


89 — 




Nd2-xCe^Cu03^ 


[35] 


21 - 




La2Cu04+5 


[36] 


32 10* - 10^ 




Ba-K-Bi-0 


[37] 


32 - 





The PME was first observed in the Bi2Sr2CaCu208-i-5 
(Bi2212) the structure of which is schematically shown in 
Fig. 2. The lattice constants are a = 5.14lA, b = 5.148A 
and c = 30.66A. The onset of superconductivity occurs 
at Tr « 87K and the transition width is about 5 K. In 



general, Tc depends on the oxygen content [4] . 
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FIG. 2. Schematic structure of the crystal lattice 
Ba2Sr2CaCu2 08+a; where triangles represent pyramids with 
an oxygen atom in each corner. The lattice constants 
a = 5.UlA,b = 5.418A and c = 30.661. Adapted from 
Ref. [39]. 



the voids of the multiply connected sample. The screen- 
ing in the ZFC mode is, however, larger than 100% of 
— l/47r. It is because the screening currents, which are 
induced by applying the field below Tc at the start of the 
measurement, screen also voids present in those granu- 
lar materials [4]. Therefore the screened volume is larger 
than the effective volume with respect to the mass of the 
sample as calculated by the x-ray density. 
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FIG. 3. ZFC and FC susceptibilities as a function of tem- 
perature of a high-quality ceramic Bi-2212 in very low fields. 

The reduced FC flux expulsion resembles the behavior com- 
monly observed at intermediate field. After Braunish et al. 
[4]. 



The Bi2212 samples are prepared by two techniques: 
the reaction sintering [2,4,38] and the melt-cast [4]. As 
shown by Scanning Electron Microscopy, sintered Bi2212 
have a porous structure with a typical grain size of 1 
- 4 jLtm [4] (see also Tab. 1) and typical thickness of 
0.2 - 0.4 /xm [39]. High resolution transmission elec- 
tron microscopy studies [33] have shown that the samples 
obtained by the melt-cast technique have an extremely 
polydomain microstructure on a /xm length scales. The 
grains are significantly larger and more densely packed 
compared to the sintered samples. The single crystalline 
domains are preferably c-axis oriented with extremely 
sharp interfaces (width < 1 nm) which provide good con- 
tacts between the domain in a6-planes. Such kind of in- 
terfaces result in Josephson junctions with quite large 
critical currents (10^ — lO^A/cm^ [33,40]). Depending on 
the morphology of the samples the Bi-2212 may exhibit 
the PME or not. 

Fig. 3 shows the standard Meissner behavior 
of the dc susceptibility in granular superconductor 
Bi2Sr2CaCu208 under FC and ZFC conditions. The 
results were obtained on a high-quality sintered sample 
[3,4]. In the FC regime one has a fractional Meissner ef- 
fect due to a certain amount of magnetic flux trapped in 




m m ' - w '""^ m. 



FIG. 4. ZFC and FC susceptibilities as a function of tem- 
perature of a granular superconductor Bi2Sr2CaCu2 08 ex- 
hibiting the PME in very low fields. Reprinted from [3]. 

The unusual behavior of the dc susceptibility in Bi- 
2212 compound obtained by the melt-cast technique is 
shown in Fig. 4. For very small external fields {H ~ 
0.01 — 1 Oe) a paramagnetic response in the FC con- 
ditions apperears below the superconducting transition 
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temperature. The smaller the H, the stronger the para- 
magnetic signal, as seen from Fig. 5. 



Ci 




-30 - 
-40it- 
-50 



91. i-.a:1:3 
* VSM 




1 



ICO 



IC 

H <0e; 

FIG. 5. Field dependence of the FC susceptibility of the 
melt-process sample KnBockl and of the sintered sample 
Ne41T10B at 20K. The applied field was varied between 10 
mOe and 10 kOe. The measurements were performed with the 
SQUID magnetometer, the vibrating sample magnetometer 
(VSM), and a Faraday balance (see [4] for details). The PME 
of KnBockl sample emerges in fields H < 2 Oe. Reprinted 
from [4] . 

For fields H > 1 Oe, x becomes slightly diamagnetic. 

The fact that even small fields arc sufficient to enforce 
nearly complete polarization reveals that orbital cur- 
rents spontaneously created in the superconducting state 
should interact weakly. It was shown that the magnetiza- 
tion consists of two parts: the conventional diamagnetic 
part and the paramagnetic one [4,41] 



M{H) 
Mo{H) 



(2.2) 



X.dia ~l~ 



(Ho + HY 



Here Xdia, Mg, a and Hq depend on the samples {Hq may 
be interpreted as the field below which the thermal and 
interaction effects suppress complete polarization). For 
Bi-2212, e.g., Xdia is about 17% of -l/47r,a ^ l,Ho = 
0.16 Oe and M,=0.095 G/47r. 

As to the ZFC experiments, at a first glance one may 
think that the Meissner effect is similar to that in ordi- 
nary ceramic superconductors. Namely, Xzfc depends 
on T but not on the external field. However, a more 
careful inspection [2,38] shows its rather weak field de- 
pendence. This nonlinear effect gets enhanced near the 
phase transition to the superconducting state. 

Recently Freitag et al. [33] have found a clear qualita- 
tive distinction between the morphology of the samples 
with and without the PME when approaching a length 
scale of several microns. It is shown in Fig. 6 where High- 



Resolution Transmission Electron Microscopy photos of 
two respective samples are shown. The samples with 
PME are characterized by an extremely polycrystalline 
" chaotic" structure with typical sizes of single crystalline 
domains < /um (see Fig. 6a). In contrast, as seen in Fig. 
6b, samples without PME show a pattern of large well- 
ordered single crystalline domains. The density of the 
domain boundaries is much larger for the former samples 
than for the latter. 





— • ifcesL " , i 

FIG. 6. Higli-RcKohitiou Transmission Electron Mi- 
croscopy photos of samples with PME (a) and without PME 
(b) on the scale of several microns. After Freitag et al. [33]. 



The another class of materials where the PME has 
been also observed is high-quality, twined single crystals 
YBa2Cu307_5 [34,42]. The study of this material would 
help to discriminate the d-wave and the grain bound- 
ary origin of 7r-j unctions. Since crystals are single the 
later mechanism can be ruled out. The PME has found 
to appear only with the magnetic field parallel to the c 
axis, i.e., the spontaneous currents are confined to the ab 
plane [34]. It is clear from Fig. 7, the magnitude of the 
observed signal here is considerably smaller than that for 
the ceramic samples (less than 3% of the full shielding 
value as compared to up to 50%) [3,4,23]. The field be- 
low which a PME is observed is, however, much higher 
than that in the ceramic samples (0.7mT compared to 
about 0.05 mT). Since the sample is a single crystal one 
has a complete flux expulsion in the ZFC regime. 



A positive FC magnetization has also been observed in 
Nd2-a;CexCuOj/ [35], La2Cu044-5 [36,43], and Ba-K-Bi- 
O compounds [37]. In the later case the magnitude of of 
Xfc is much smaller than the one observed in BSCCO 
samples at weak fields but it remains positive for much 
larger applied fields. 



Table 1 collects ceramic high- Tc superconductors show- 
ing the PME. 
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B. Conventional superconductors 

One of the scenarios for the PME in granular hole- 
doped cuprates discussed above is that the FC paramag- 
netic response appears from current loops with tt phase 
shifts of the superconducting order parameter at some 
grain-boundary junctions. It was argued that such be- 
havior would be expected to occur in a d-wavc super- 
conductor, but not in a conventional s-wave supercon- 
ductor. The test of this hypothesis led to the discov- 
ery of the PME in conventional superconductors like Nb 
[44,16,17,45], Al [18] and the Nb-AlO^-Nb tunnel junc- 
tions [46,47]. The occurrence of the PME in these sys- 
tems is probably due to the confined geometry and flux 
trapping. In this section we discuss the main experimen- 
tal findings for the conventional superconductors. 



0.0 

-0.4 
-a .6 

-o.a 

-1.0 



- (a) 



0.02 n-T 
0.? mT 
0.5 mT 
1 .0 mT 



m 
iti 

? ■ « 



t 

nnfYiTiiliiriTiril'---'' 



s 



O.OJ 
0.Q2 
0.01 

o.oa 

■0.0^ 



SP 90 91 


92 


T (Kl 




. 1 - ■ 1 






0.02 mT 




0.2 mT . 




0.5 mT 


\ n 


1 .0 mT 










O Q C ". ^ O C L-- C 




1 . 1 . . 





£3 



fl9 



92 



M PI 
T |Kf 

FIG. 7. Susceptibility of single crystal YBa2Cu307-5 in 
various applied fields vs temperature T. Prom Riedling et al. 
[34]. 

In the experiments of Minhaj et al. [44,16] Niobium 
disks of diameter 6.6 mm and thickness 0.127 mm were 
used to study the PME. The dc FC magnetization shows 



the paramagnetic response below the superconducting 
transition temperature Tc ~ 9, 2 K when the applied 
magnetic field is normal to the disk. It is demonstrated 
in Fig. 8. The paramagnetic signal is much weaker com- 
pared to its counterpart in BiSrCaCuO compounds and 
it becomes diamagnetic for fields larger than 20 Oe. 
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FIG. 8. The FC (upper) and ZFC (lower) susceptibility for 
a macroscopic Nb disk with magnetic fields applied normal to 
the disk surface. The data are scaled to the complete shield- 
ing value of xzFC. The inset shows an enlargement of xzfc 
around T = 9, 26 K. Prom Thompson et al. [16] 

Comparing Figs. 7 and 8 one can see the remark- 
able similarity in the temperature dependence of the sus- 
ceptibility of Nb disks and single crystal YBa2Cu307_5. 
With increasing temperatures the ZFC data indicate the 
presence of two different superconducting transitions at 
Ti w 9.2K and Tp « 9.06K. The strong field dependence 
of Xzfc indicates that the local field is larger than the 
lower critical field Hci{T) in the interval Tp < T < Ti. 
The appearance of the two transition temperatures is 
also seen in the FC measurements: at Ti the paramag- 
netic moment first appears (vanishes) and a lower tem- 
perature Tp defines the temperature where the positive 
moment no longer increases [16,48]. The sharp increase 
of positive FC magnetization upon cooling below T\ is 
fairly spontaneous similar to the onset of global diamag- 
netic screening currents at T^ rather than the viscous 
nature exemplified by flux flow. Below Tp this moment 
apparently does not change with temperature and can 
be regarded as an additive constant to the FC magnetic 
moment of a non-PME superconductor. 

To gain more insight on the nature of the FC paramag- 
netic moment for the Nb disks, magnetic hysteresis loops 
were measured [16,45]. A set of such loops recordered in 
the temperature range from 9.04 to 9.09 K, i.e., around 
Tp are shown in Fig. 9. Between Tc and Ti the system 
is superconducting but only a small diamagnetic screen- 
ing current can be induced by ramping the external field. 
Below T\ the critical current is expected to increases sig- 
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nificantly and the loops exhibit a strange, nearly paral- 
lelogramlike shape the size of which increases with lower- 
ing temperature [45]. Just below Tp the hysteresis loops 
change their shapes drastically and they arc more rem- 
iniscent of a type-II superconductor with the magneti- 
zation becoming less diamagnetic as the flux penetrates 
into the bulk for H > Hci- 
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FIG. 9. Magnetic hysteresis loops of the Nb disks in the 
temperature range 9.04 < T < 9, 09 K for 0.01 K increments. 
The remarkable change of their shapes between 9.05 and 9.06 
K, i.e., at T = Tp. Prom Pust et al. [45]. 

Interesting similarity to the YBa2Cu307_a; single crys- 
tals [34] results in the hysteresis behavior of the FC mag- 
netization between cooling and warming cycles. As one 
can see from Fig. 8 the FC magnetizations are more 
positive than its warming coimter part. On cooling, the 
flux should overcome activation-type process in order to 
be expelled and consequently the FC magnetization is 
more positive than the equilibrium one. In the warm- 
ing cycle, the flux cannot, however, easily penetrate the 
sample and the magnetization is more diamagnetic. This 
behavior is common in type-II superconductors and is 
consistent with a theoretical interpretation based on the 
critical-state model [49]. 

It has been found that the PME in Nb samples is very 
sensitive to their surface and geometry [16,17,45]. In par- 
ticular, the paramagnetic response occurs in Nb disks 
cut from some sheets of rolled Nb. Disks cut from other 
source materials do not show the PME [17]. Fig. 10 
shows the FC magnetizations obtained before and after 
polishing the sample faces. Interestingly, the FC param- 
agnetic signal disappears after both faces of the sample 
were polished. This gives strong evidence that the sur- 
face pinning of the magnetic flux plays an important role. 
One may think that the PME is likely to arise from in- 



homogeneously trapped flux, and is unlikely to have any 
relationship with d-wave superconductivity. 
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FIG. 10. FC (a) and ZFC (b) magnetizations of a Nb disk 
measured before (solid squares) and after (solid triangles) pol- 
ishing the sample surface. The values of the external magnetic 
field are shown next to the curves. Inset: a high point density 
ZPC measurement at 0.05 Oe from another sample cut. The 
PME disappears after both faces of the sample were polished. 
Prom Kostic et al. [17]. 

We now disciiss the PME in small (micromctrc-size) 
superconductors. The experiments were carried out for 
Nb and Al discs [18] with diameters from 0.3 to 3 /im and 
thicknesses from 0.03 to 0.15 fim. The thickness is com- 
parable to the characteristic penetration depth A ~ 500A 
[50]. Study of small samples is interesting in two ways. 
First, an application of the Koshelev-Larkin model [19] to 
data for Nb discs [17] suggested that the paramagnetic 
magnetization must be small for samples much thicker 
than A; however, for very thin films one can expect 
a much larger effect due to macroscopic penetration of 
Meissner currents into the sample interior. Second, con- 
finement of superconductivity in a small volume leads to 
pronounced quantization , so that a mesoscopic super- 
conductor resides in one of well-resolved states, depend- 
ing on temperature and magnetic field (superconducting 
states are characterized by a different number and distri- 
bution of vortices [51-54]). This allows one to measure 
magnetization of individual vortex states and rules out 
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the pinning effect. 
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FIG. 11. The field dependence of the Meissner response for 
Al disks of diameter 1.0 (open squares) and 2.5 /umi (closed 
circles) at T=0.4 K. The dashed line is a guide to the eye. 
The inset compares FC and ZFC magnetization of the 2.5 
/um disk at the field where the paramagnetic response is close 
to its maximum. The jumps in the ZFC curve correspond to 
entry of individual vortices into the disk interior. After Geim 
et al. [18]. 



The typical field dependence of the FC paramagnetic 
moment of mesoscopic samples is shown in Fig. 11 [18]. 
The strong oscillation clearly seen for the larger sample 
of diameter 2.5 /xm is due to size quantization. Each 
jump corresponds to a change in the number of vor- 
tices inside the disk, which can either form an array 
of single quantum vortices or assemble into a single gi- 
ant vortex [21,51 54]. The latter configuration is ex- 
pected at fields between the second and third critical 
field, Hc2 < H < Hcs, where Hcs is related to the sur- 
face superconductivity [21]. The smaller sample (diam- 
eter 1.0 /im) docs not show the paramagnetic response 
over the entire field interval. This qualitatively different 
behavior may be explained by the fact that, in smaller 
sample the superconductivity is suppressed by ~ 3 flux 
quanta entering the whole disk area while ~ 20$o ($o 
is the flux quantum) are needed to spoil superconductiv- 
ity of the larger disk [18,51-53]. Thus, the diamagnetic 
response is always observed in low magnetic flelds and 
the PME occurs only in intermediate flelds allowing at 
least several flux quanta penetrate into the interior. This 
seems to contradict the previous studies on macroscopic 
samples, where the paramagnetic signal is recorded in 
very low flelds and gradually disappeared with increas- 
ing field. Such discrepancy is due to the large number of 
flux quanta (many thousands) may enter the macroscopic 
samples even in lowest flelds. 
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FIG. 12. The magnetization measured by cooling in a field 
and by sweeping the field at a constant temperature (0.4 K). 
Arrow shows the direction of the sweep. The FC data shown 
by open squares are for the 2.5 ^m disk of Fig. 11. The 
solid curves were measured by pausing at various fields and 
then sweeping the field up and down. When the field swept 
continuously, the magnetization evolves along one of the solid 
curves until it reaches the end of this curve and jumps to the 
next one belonging to another vortex state. Then, the pro- 
cess repeats itself all over again as shown by dashed lines. The 
closed squares at 35 G indicate the low-temperature states for 
the FC and ZFC curves of Fig. 11. The inset illustrates the 
compression of a giant vortex into a smaller volume which al- 
lows extra flux to enter the sample at the surface. From Geim 
et al. [18]. 

The nature of the PME becomes more evident from 
measurements of the magnetization by sweeping the mag- 
netic field at a constant temperature [18] (sec Fig. 12). 
Instead of a single magnetization curve characteristic for 
macroscopic samples one has a family of curves corre- 
sponding to different vortex states. Since several super- 
conducting states are possible at the same magnetic field, 
the states in which the PME occurs are metastable. Only 
the state with the most negative magnetization is ther- 
modynamically stable [52,53]. The metastability results 
from the sample surface which favours a large supercon- 
ductivity and opposes both vortex escape and entrance 
[55] . As shown below, combining the ideas of flux com- 
pression [19] and the trap in giant vortex states [21] one 
can explain the experimental findings for mesoscopic su- 
perconductors [18]. 



III. MECHANISMS FOR n JUNCTION 

A. Scattering on paramagnetic impurities 

The effect of magnetic impurities present in the barrier 
between s-wave superconductors on the tunneling cur- 
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rent was first studied by Kulik [56] (sec also [57]). Later 
Bulaevskii et al. [58] have reahzed that the presence of 
such impurities may lead to a spin-flip intermediate state 
which actually shifts the phase by tt. 

Consider a SIS contact with magnetic impurities in the 
thin insulator layer. The Hamiltonian of this system is 
as follows [56,58] 

+ c.c. (3.1) 

Here cij:^{bj:^) are electron annihilation operators in su- 
perconducting layer A (B); Sft is spin of impurities local- 
ized at point n and a are Pauli matrices. 

Using the Green function technique [59] , in the second 
order of perturbation theory we obtain the stationary 
Josephson current J and the contact energy E in the 
following form 

h 

J = Jc sin (f), E — — — Jc cos (j), Jc = Jo — Js 
2e 

Jo = 2w^et^N^{0)A{T)th^, 

J, = 27r'eY,vlS{S + l)N\0)A{T)th^, (3.2) 

where t"^ and are average values of i-- and ^ on the 

kk' kk'n 

Fermi surface respectively, A^(0) - the density of states 
on the Fermi surface, 2A(T) - the energy gap and - 
the difference between superconducting phases. Jg is the 
standard Ambegaokar-Baratoff term and Jg describes the 
electron tunneling with the spin flip [56] . 

Since ugg,^ ~ [ed(ed (7)]"^ [60], where td - energy of 
one electron and U - the energy of the Coulomb repul- 
sion between two electrons, one can choose so that Jc 
becomes negative (Js > Jo) [58]. The another possibility 
to have Jc < is that, as follows from Eq. (3.2), the con- 
centration of magnetic impurities should be high enough 
making the spin flip term to be dominant. 

Rewriting the Josephson current as 

J _ j Jc sm{(f) + 0) if Jc > , . 

\ |Jc|sin(0 + 7r) if .J, < ^ ' 

one can understand the meaning of terminology tt- 
junction. Namely, in the case when the critical current 
Jc < the current (or the contact energy) has the same 
form as in the conventional case with Jc > [10] but 
the phase is shifted by tt. For this reason an junction 
with negative critical current is called tt junction whereas 
Jc > corresponds to the standard junction. 

An another possibility of occurrence of negative criti- 
cal currents was also pointed out by Spivak and Kivelson 
[61]. The main difference from the Bulaevskii approach 
[58] is that they described an impurity by the Anderson 



model in a regime where the localized magnetic moments 
exist. Spivak and Kivelson argued that Jc may become 
negative in a dirty system due to an interplay of disorder 
and electron correlations. Similar results were obtained 
earlier by Altshuler et al [62]. 

The existence of 7r-junctions leads to non-trivial physi- 
cal phenomena. In the next section we will show that the 
spontaneous current in the ground state of a frustrated 
loop of odd number of tt junctions is nonzero. This can 
serve as a reason for the paramagnetic response to the 
external magnetic field. 



B. d-wave picture of a 7r-junction 

In all cases discussed above the appearance of the tt 
junction depends not on the intrinsic properties of a su- 
perconductor, but rather on external conditions like de- 
gree of disorder, impurity concentration etc. Further- 
more, the spin flip mechanism, for example, may re- 
quire the high density of magnetic impurities causing the 
strong magnetic interaction which would invalidate elas- 
tic tunneling. 

There is, however, another, intrinsic and more exciting 
scenario for Tr-junctions: the spontaneous moments may 
naturally occur in ceramic superconductors with the d- 
wave pairing symmetry [9]. This idea was inspired by 
the the works of Geshkenbein and co-workers [63,64] on 
heavy-fermion superconductivity showing that special ar- 
rangements of junctions may lead to the tt shift. In 
this chapter, using the Ginzburg-Landay theory for the 
Josephson jimction with two connected li-wave supercon- 
ductors we show that the critical current may be positive 
or negative. 

It is well established that the phonon mediated 
electron-electron interaction gives rise to spin-singlet 
pairing with s-wave symmetry [65] in conventional low- 
Tc superconductors. The question of symmetry of the 
order parameter of high- Tc superconductors is under in- 
tense debate [66,67]. Although the mechanism of high- 
Tc superconductivity remains unknown experiments on 
Shapiro steps [68], the magnetic-flux states of YBCO 
[69], the Andreev-reflection [70] and spin susceptibilities 
[71,72] give strong evidence for spin-singlet Cooper pair- 
ing in ceramic materials. We shall focus on this type of 
pairing. 

Within the Ginzburg-Landau formalism the normal 
metal - superconductor phase transition may be de- 
scribed by the scalar complex order parameter ^'(r) 

which is the pair wave function. In momentum-space 
representation, \l/(fc) oc< aj:^a_j^^ > and is related to the 

gap function through ^'(fc) = A^/2E^, where Ej: is 
the quasiparticle excitation energy. The energy gap A^ 
can serve as a well defined superconducting order pa- 
rameter or as a measure of long-range phase coherence 
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in the pair state. Its symmetry can be experimentally 
determined, even without detailed knowledge about the 
microscopic origin of superconductivity [15]. 

Since the transition to the superconducting state is sec- 
ond order it should be accompanied with a continuous 
symmetry breaking. The symmetry group H describing 
the pair state must be a subgroup of the group Q describ- 
ing the normal state, i.e. H C G [73]. Considering pairing 
in a crystal one has Q = Gx Rx Ix Tx U(l), where G is 
the finite crystallographic point group, R the symmetry 
group of spin rotation, I the space inversion operation, 
T the time-reversal symmetry operation and U(l) the 
one-dimcnsional global gauge symmetry. For BCS su- 
perconductors the gauge symmetry U(l) is broken in the 
superconducting state due to the existence of off-diagonal 
long-range order. In an unconventional superconductor, 
in addition to U(l) one or more symmetries may be bro- 
ken at Tc- Since spin-orbit interaction in ceramic super- 
conductors is rather weak the pair wave function should 
be either singlet (spin 5 = 0, fc) = \E'(A;)) or triplet 
{S = 1, = The mixed state like a^^+b^* 

is forbidden. Due to the spatial inversion the singlet pair 
wave function does not change under the inversion op- 
eration {I^^ = whereas the triplet one changes its 
sign = -vp*). 

Point group classification of pair states in cuprate su- 
perconductors has been extensively studied [66,74-79]. 
They arc divided into two groups: tetragonal crys- 
tal lattice with point-group symmetry D4fi and or- 
thorhombic crystal lattice with point-group symmetry 
D2h- La2_:rSr^Cu04, TlaBaaCaCuaOs, HgBa2CaCu04 
and some YBa2Cu3 compounds belong to the first group 
whereas YBa2Cu307_5 and Bi2Sr2CaCu208 possess D2h 
symmetry. All possible pair wave functions of D^h and 
D2h groups are given in [76]. It should be noted that 
superconductivity of cuprates basically originates from 
the Cu02 square/rectangular layers. Therefore we focus 
on the symmetry of possible pair states on the square 
lattice which is characterized by point-group symmetry 
Ciy. The list of spin-singlet even parity pair states for the 
tetragonal point group Ciy is given in Tab. 2. All repre- 
sentations Aig,A2g,B-ig and B2g are one-dimensional. 

TABLE 2. Even-parity pair states allowed by sym- 
metry 



Name State gap function A(fc) Nodes 



s-wave 


'Alg 


1 


none 


9 


\A2g 


kxky{k'^ — ky) 


line 




^Blg 


line 


^xy 


^B2g 


kxky 


line 



Since there is strong experimental evidence for dx^-y^ 
symmetry [80,15] in the high-temperature superconduc- 
tors (see below) we will concentrate our study on this 



type of symmetry. Neglecting the possible orthorhom- 
bic distortion the pair wave function is described by Big 
representation of or D^h group. In momentum space 
the pair wave function is 

^{k) = < cg^c_gj^ > = coskx — cos ky . (3.4) 

The schematic description of the pair wave function 

of s and dx2^y2 symmetry is shown in Fig. (13). The 
anisotropic function cos kx — cos ky has the same symme- 
try property in C^y as k'^ — ky. For this reason, this state 
is called the d-wave with the angular momentum L = 2, 
although the classification with respect to angular mo- 
mentum has no real meaning under the discrete crystal 
field symmetry. Since function (3.4) is described by one- 
dimensional representation, the time-reversal symmetry 
can not be broken at Tc. Breaking of this symmetry is 
possible only for multi-dimensional representations [81]. 

We now consider the phenomenological theory of the 
Josephson effect in the case of a junction between two 
unconventional superconductors based on symmetry ar- 
guments. So far as dx2-y2-^Na\e corresponds to one- 
dimensional representation B\g their complex order pa- 
rameters r?i and r/2 should be one-component 

r?i = |r?i|e'^Sr?2 = |r?2|e''^^ . (3.5) 

An intensive discussion of the experimental evidence for 

XY-universality class identifying the order parameter as 
a complex scalar belonging to a one dimensional irre- 
ducible representation may be found in Ref. [82] . 




d-wave 



FIG. 13. Schematic description of s- and dj,2_j^2-wave func- 
tions. The order parameter may be positive (-I-) or negative 
(-)■ 

The free energy of two subsystems (see Fig. 14) joined 
by a junction is determined by the expression [10,83] 

= ^1+^2+ T\2 , 

Ti = { drFi (i = l,2), 

Jrni>0 

KailDxVil'^ + \Dyrjif) + Ki2\D,r]if + — , 
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T\i = j> dS togi{ni)g2{n2){il*V2 + ?7i»?2) • (3-6) 

Here ao, P, Kij(i,j = 1,2) and to are phenomenologi- 
cal parameters; outward normal to the junction surface 
vectors fii [fii = —712) arc defined in the basis of local 
crystal axes; D^^y^z are the components of the gauge- 
invariant gradient 5 = (V — 2mA/ where A is the 
vector potential (B = V x A). Ti and T2 arc scalar un- 
der the operations of the complete symmetry of the sys- 
tem Q including the crystal symmetry D^h (or C^v), time 
reversal, spin rotation and gauge symmetry [84,75,76]. 
The interface term Ti2 was chosen in the lowest order 

5A 

-r ^^2, or 

SA 



and so that it generates the correct boundary conditions 
[63,75,85,84]. Functions (?i(ni) and 52(^2) should be in- 
variant under all operations of the symmetry group of 
the respective half-space. This condition is satisfied if we 
choose them to have the same symmetry properties as 
the corresponding order parameters. For the d^^-yz we 
have 

gi{fii) = rii^ - riiy ~ cosnj^ - cosn^y {i = 1, 2). (3.7) 

Varying the gradient term of Eq. (3.6) with respect 
to A, we obtain the current density in the left and right 
half-spaces 



ad 



2ttc 

-—Rer]t[Kn{D^ + Dy) + Jfia-Dz]??! , 
27rc 

-—Rev*2[K2i{D^ + Dy) + K22£'.]?72 • 
$0 



(3.^ 



The boundary conditions for the Ginzburg-Landau ing the full free energy (3.6) with respect to r/* and 772 
equations at the junction interface arc obtained by vary- 

i[Kii{nixD^^ + niyDy) + Ki2nizDz]r]i\s = -io5i("i)32(«2)??2|s , 

i[K2i{n2xD:c + n2yDy) + K22n2zDz]r]2\s = -to9i{ni)g2{n2)r]i\s ■ (3.9) 



From Eqs. (3.8) and (3.9) we have the current density 
though the Josephson junction 

J = (Ji+j2)ni|s = Jc sin(0i - (/)2) , 

= ^!^g-i^(fii)g2{n2)\m\\m\ ■ (3-10) 

Taking Eqs. (3.7) and (3.10) into account one can see 
that the critical current Jc oc (nf^, — n\^){n{y — n\y) may 
be positive or negative depending on orientation of vec- 
tors ni and 712 . A negative J,, is equivalent to an intrinsic 
phase shift of tt in the jmiction. Because gi{ni) has es- 
sentially the same symmetry as the order parameters, the 
Josephson effect allows one to probe the phase of the pair 
wave function. The existence of tt junctions may serve as 
an indicator for the d-wave superconductivity. 




FIG. 14. Josephson junction between two rfj,2_j,2-wave su- 
perconductors. Vectors ni and n2 defined in the basis of local 
crystal axes are normal to the interface. 

C. Experimental verification of d-wave pairing 

So far as the c?-wave pairing may play a key role in 
understanding the PME in ceramic superconductors, in 
this part we discuss experiments which support its exis- 
tence. Those experiments are divided into two classes: 
phase insensitive and phase sensitive ones [15]. 
Phase insensitive experiments. There are two kinds of 
phase-insensitive experiments. One of them probes the 
gap features at the Fermi surface. It includes, in par- 
ticular, the angle- resolved [86], tunnel [87], Raman scat- 
tering [88] and point-contact [89] spectroscopies. The re- 
sults obtained by these techniques for a number of specific 
hole-doped oxides, are usually interpreted as manifesta- 
tion of the d-wave pairing although it is impossible in 
some cases to distinguish it from the anisotropic s-wave 
symmetry [15]. Here we focus on the s(;cond group of 
phase-insensitive experiments based on the fact that the 
low temperature behavior of thermodynamic and trans- 
port properties of d-wave superconductors is governed 
not by exponential as in the BCS theory but by power 
laws. 

At low temperatures the specific heat of of BCS super- 
conductors behaves as Cs{T ^ 0) ~ T'^'^ exp(-Ao/T), 
where Aq is the isotropic s-gap at T = [7]. Simple 
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calculations [90] show, however, that for the hexagonal 
dx2-y2 superconductors ~ due to lines of nodes on 
the Fermi surface. There is some controversy about ex- 
periments of different groups. The Cs ^ was found in 
Bi-compounds [91] whereas the linear behavior Cs ^ T 
was reported by other investigators [92-94]. Recent ex- 
periments on single crystals [95,96] and ceramic samples 
[97] of YBCO show, however, that ~ T^. 

In the presence of the magnetic field the quadratic tem- 
perature dependence rolls over to a T [98], i.e. C,, = 
k-fnT^H/H^2 for T/Tc < ^/W/H^, where is a con- 
stant of order one and 7„ is the coefficient of the linear-T 
term in normal state [7]. This linear-T and H'^/'^ depen- 
dence was confirmed by experiments [96,97,99]. Further- 
more Revaz et al [100] found that the anisotropic com- 
ponent of the ficld-dcpcndcnt specific heat Cs{T,H\\c) — 
Cs [T, _L c) of single crystal YBCO obeyed a scaling re- 
lation predicted for a superconductor with a line of nodes 
[101,102]. 

In the low temperature limit the London penetration 
depth Xl{T) of the conventional superconductors ap- 
proaches its zero value Al(0) exponentially: SXl{T) = 
(Al(T) - Al(0))/Al(0) - T-i/2exp(-Ao/r). The the- 
ory of dx2-y2 superconductivity gives 6Xl{T) T for 
pure materials and S\l(T) ^ for dirty ones [103,104]. 
For nominally clean YBCO samples linear dependence 
was observed [105] whereas for Zn- and Ni-doped as well 
as nonhomogeneous crystals the quadratic behavior was 
reported [105-111]. Although these results agree with 
the d-wave symmetry of the gap function, they cannot 
distinguish between the d- and s-wave because the linear 
dependence of JAl(T) could also arise from promixity 
effects between alternating s-wave superconducting and 
normal layers in the cuprates [112,113]. 

The another way to probe the symmetry of the or- 
der parameter is to study the low-T behavior of the 
thermal conductivity of electrons Kp{T). Ke{T) ^ 
exp^—Ao/T) for the BCS superconductors but 
KeiT) ~ T for the d^2_y2 materials [114,115]. This 
linear dependence was shown to be insensitive to the 
impurity scattering [115]. Unfortunately, it is difficult 
to extract the thermal conductivity component Ke{T) 
from experiments due to the complex mutual action of 
electrons, phonons and impurities [114,116]. Neverthe- 
less, the existing experiments confirm the linear temper- 
ature dependence in Zn-dopcd YBCO [117-119], below 
T* = 200mK in Bi2Sr2Ca(Cui_a;Ni:r)208 [120] and in 
BiSrCaCuO [121]. 

As to the ultrasonic attenuation coefficient ag in su- 
perconducting state with anisotropic pairing, the the- 
ory [122] predicted ag ~ T at low temperatures. Ex- 
periments performed for YBa2Cu307_x [123,124] and 
Lai.8Sro.2Cu04_a; showed that as ~ T" with a large 
scatter of exponent n for each substance. 

Finally, as T ^ 0, the nuclear spin- lattice relaxation 
rate Tj~^ of the dx2-y2 superconductors was shown to 



scale with temperature T^^ ^ [104]. Such a 
temperature dependence was confirmed by experiments 
[104,125] 

It should be noted that low temperature asymptotics 
of various thermodynamic and transport properties of s- 
wave superconductors can roll over power-laws [126,127] 
provided the gap values are widely distributed due to 
structure domains, charge stripes, charge-density waves 
or other mesoscopic inhomogeneities. So the interpreta- 
tion of the rf-wave symmetry based only on the low-T 
asymptotics remains ambiguous. 

Phase sensitive experiments. Recently, a number of 
phase-sensitive experimental techniques has been devel- 
oped to determine the order parameter symmetry of 
cuprate superconductors [80,104,66,15] starting from the 
pioneering work of Wollman et al [128]. The key idea is 
relying on the sign changes in the Josephson critical cur- 
rent Ic as proposed first by Geshkenbein and Larkin [63] 
as a test for axial p-wave pairing symmetry in the heavy 
fermion superconductors. 

In the "corner SQUID" design shown in Fig. (15a) 
Josephson junctions are made between s-wave Pb thin 
films and two orthogonally oriented ac and be plane faces 
of single crystal YBCO [128]. If YBCO is a d-wave su- 
perconductor, then, according to Eq. (3.10), there should 
be a TT phase shift between weak links on adjacent faces 
of the crystal. For the "edge SQUID" geometry in Fig. 
(15b) no phase shift is expected because both two junc- 
tions are on the same crystal face. 

We consider the dc "corner" SQUID with junction crit- 
ical currents la and J^. Then a bias current through the 
SQUID is the sum of the currents passed by junctions a 
and h 

I = la sin <j>a + h sin (pb . (3-11) 

On the other hand, the quantum phase coherence around 
the SQUID loop forces the constraint on the gauge- 
invariant phase differences across the junctions on the 
ca and cb faces of crystal [7,128], namely, 

</)„-0b + 27r^^-^ + ^J+5„b = 27rn, 

(3.12) 

where La and Lt are the effective self-inductances of the 
two arms of the ring, n is integer, 6ab accounts for the in- 
trinsic phase shift inside the YBCO crystal between pairs 
tunneling into the crystal in the a and c directions. One 
can expect that Sab = ti" and 5ab = for the " corner" and 
"edge SQUID", respectively. 
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FIG. 15. Experimental geometry for the SQUID and single 
junction experiments (see Ref. [128,137]): (a) Configuration 
of the corner SQUID interferometer experiment used to deter- 
mine the relative phase between orthogonal crystal directions 
a and b. (b) The edge SQUID used as a control sample, in 
which both junctions are on the same crystal face, (c) The 
corner configuration for a single junction experiment, (d) The 
same as in (c) but for the edge geometry. 

In order to find the maximum current one can substi- 
tute (pb from Eq. (3.12) into Eq. (3.11) and then find the 
maximum of / with respect to 4>a- For a symmetric dc 
SQUID with equal junction critical currents la = lb = lo, 
in the limit of zero loop inductance, simple algebra gives 
the following expression for the field dependence of the 
maximum supercurrent 



ext) 



= 2/n 



cos TT 



Sab 

2 



(3.13) 



For finite values of self-inductances the modulation 
depth of the critical current is reduced because the cir- 
culating currents generate a flux contribution in the ring 
[128]. From the last equation it is clear that if YBCO 
has s-wave symmetry then 5ab = and the circuit will 
behave as an ordinary dc SQUID [7] : Ic has a maximum 
at ^ext — 0. Such a behavior should be also exhibited 
by "edge SQUID" regardless to the type of symmetry of 
pairing in the YBCO material. In contrast, for rf-wave 
symmetry the "corner SQUID" should have 6ab = tt and 
the critical supercurrent would display a minimum at 
zero external flux. This important prediction was con- 
firmed by the experiments of Wollman et al [128]. 

It should be noted that there are several complicat- 
ing factors in interpretation of results obtained by the 
SQUID experiment of Wollman et al [80,15]. Since 
YBCO has an orthorhombic crystal structure, it has a 
tendency to form a twin boundaries at which a and b lat- 
tice constants are interchanged. This would randomlize 



the phase that these experiments depend on. However, 
there are both experimental [129,130] and theoretical [80] 
showed that the order parameter maintains its orienta- 
tion across twin boundaries and forms a single domain 
even in twined samples. Moreover, experiments on de- 
twined samples [131,80] gave consistent results indicat- 
ing that the d^2_y2 component of the order parameter 
has the same phase across twin boundary. 



Klcmm concerned [132] about the role of the corners 
in dc SQUID measurements of the phase anisotropy of 
YBCO crystals [128]. Since flux trapping, demagnetiza- 
tion and field-focusing effects can be strongly dependent 
on the sample geometry he argued that the n phase shifts 
seen between the corner and edge SQUIDS could result 
simply from their geometry differences, even for s-wavc 
superconductors. This point of view was criticized by 
Wollman et al [133] showing that corners play no signifi- 
cant role. 



An another concern is the flux trapping near the 
SQUID. This flux can be coupled to the SQUID loop and 
create a shift in the flux modulation pattern that is in- 
distinguishable from the intrinsic phase shift. In order to 
avoid this effect one can cool the SQUID many times to 
get the lowest-energy state, which will be one with no flux 
trap. Magnetic imaging of the ac and be plane faces of 
cuprate superconductors [134,135] shows that there can 
be vortices trapped between the planes, even when they 
are cooled in a very small fleld. Such trapped vortices 
could affect the measured critical-current vs applied fleld 
characteristics of the SQUID's but not substantially. 



Since Eq. (3.13) was obtained for a symmetric (/„ = 
lb) SQUID with zero sclf-inductancc. Asymmetries be- 
tween critical current in the junctions and the non-zero 
loop inductance can modify the critical current pattern 
significantly. Of particular importance here is a shift in 
the pattern, since this could mimic or obscure phase shifts 
arising from the pairing symmetry. Since the asymmetry- 
induced shifts of the Id^ext) pattern are proportional to 
the screening, it is vital to design SQUIDs with negli- 
gible inductance. Such a problem was solved by Schulz 
et al [136] using YBCO thin films epitaxially grown on 
bicrystal and tetracrystal substrates. They obtained, as 
one can see from Fig. (16), nearly ideal dependences of 
the critical current on applied fleld, with a minimum at 
i? = for the 7r-ring SQUID, as expected for a d-wave su- 
perconductor in the tetracrystal geometry used. A more 
detailed discussion on the SQUID technique was given in 
Ref. [15]. 
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FIG. 16. Field dependence of the critical current of the 
TT-SQUID (a) and (b) and of the standard SQUID (c) and 
(d) at T = 77. The dip at zero field in the corner SQUID 
is evidence for d^2_y2 pairing symmetry in the YBCO. After 
Schulz at al. [136]. 

We now discuss single-junction modulation experi- 
ments [137] the geometry of which is shown in Fig. 15c 
and 15d for the corner and the edge case, respectively. 
These experiments have a distinct advantage over the 
SQUID ones since they are less sensitive to flux trap- 
ping and sample asymmetry. The basic idea is relied 
on the well known effect of an applied magnetic field on 
the critical current of Joscphson tunnel junctions [7,138]. 
The magnetic field penetrating though the barrier region 
transverse to the tunneling direction forces a gradient 
in the phase of the order parameter across the width of 
the junction, resulting in a variation of the local current 
density and a reduction in the total critical current. If 
the junction is rectangular and homogeneous with the 
field applied parallel to one edge of the rectangle, and 
if the junction size is much smaller than the Josephson 
penetration depth, the critical current has the standard 
Praunhofer form [7,138] 



In 



sin(7r<I>/$o) 



(3.14) 



7r$/$o 

which is familiar from single-lit optical diffraction, $ is 
the flux threading the junction. The Fraunhofer pattern 
given by the last equation is valid for an edge junction 
with arbitrary pairing symmetry and for a corner junc- 
tion with xisd-wave superconductors. Similar to the edge 
SQUID, this pattern has a maximum at $ = 0. 

In the case of the corner junction, the order parameter 
in the a and b directions would be of opposite sign, mod- 
ifying the diffraction pattern. In a symmetric junction 
with equal geometries on the a and b faces, the critical 
current modulates according to [137,139] 



sin2(7r$/2$o) 



7r$/2$o 



(3.15) 



Contrary to the edge junction, 7c vanishes at zero applied 
field as the current through two orthogonal faces cancels 
exactly. 
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FIG. 17. Field dependence of the critical current of the 

edge junction (a) and of the corner junction (b). The min- 
imum at zero field in the corner edge is evidence for d^2_y2 
pairing symmetry in the YBCO. After Wollman at al. [137]. 

The field dependence the critical current of the edge 
and corner junctions obtained by Wollman et al [137] is 
shown in Fig. (17). Similar patterns were reported by 
different groups [139-141]. The pattern shown in Fig. 
(17b) for d-wave corner junction does not agree closely 
with the ideal expression (3.15). It may be attributed, 
as shown by theoretical predictions, to the flux trap in 
the sample [142] or by the size of the junction [143]. So 
the experiments on the modulation pattern of the criti- 
cal current in the single Josephson junction support the 
existence of the rf-wave pairing symmetry. 

To summarize, the phase sensitive, along with phase 
insensitive experiments give strong evidence in favour of 
predominantly d-wave pairing symmetry in a number of 
cuprates. The identification of this type of symmetry is 
based purely on very general principles of group theory 
and the macroscopic quantum coherence phenomena of 
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pair tunneling and flux quantization. It docs not, there- 
fore, necessarily specify a mechanism for high tempera- 
ture superconductivity. 

IV. PARAMAGNETIC MEISSNER EFFECT IN 
D-WAVE SUPERCONDUCTORS 

This section is focused on the d-wave mechanism of the 
PME in the framework of the single-loop and multi-loop 
models. 

A. One-loop model 

As shown above, the phase across a Tr-junction between 
two d-wave superconductors is shifted by tt and the cor- 
responding energy is proportional to cos(0i — 4>2 + i^)- 
The system tries to minimize its energy by setting the 
phase difference {(f)\ — 4>2) equal to n. Thus, there is no 
way to measure this phase shift directly because it merely 
corresponds to the phase change in one of the two super- 
conductors, say, (pi ^ <pi + TT. This transformation is 
equivalent, as one can see from Eq. (3.7), to the ex- 
change of the X and y coordinates in superconductor (1). 
In other words, the connection of two d-wave supercon- 
ductors by a single 7r-junction does not, by itself, lead to 
any special observable effects. In this sense, whether a 
junction is 0- or 7r-junction is just a matter of convention. 




(a) (b) 



FIG. 18. (a) Schematic plot of a contact of three grains 
of superconductors with unconventional pairing with differ- 
ent orientation of their local crystal axis. These grains 
form a frustrated loop with nonzero spontaneous current 
in the ground state, (b) Single 7r-junction loop with the 
self-inductance L, area S and critical current 7c. 

Physically interesting consequences are, however, ex- 
pected to arise if an odd number of 7r-j unctions are con- 
nected in a loop [10]. A cartoon of a contact from three 
junctions is shown in Fig. 18a. One can choose the 
coordinates x and y in two segments to convert two tt- 
junctions into two 0-junctions but there is no transforma- 
tion to remove the remaining 7r-junction without putting 
one of the two 0-junctions back to a 7r-junction. After all 



redefinitions of the crystal axis the multi-connected loop 
has at least one 7r-junction. Since there is no way to min- 
imize the energy of all junctions and at the same time to 
keep the phase of the order parameter constant in each 
segment, such a loop is, in analogy with spin glass physics 
[144], frustrated . In what follows any single loop is sup- 
posed to contain effectively only one tt junction as shown 
in Fig. 18b or one junction which is not frustrated. 

Assuming that the loop has self-inductance L and the 
current / flowing in it is small compared to the critical 
current Ic> 0, the energy is given by the following form 
[9,10] 

i^(^,^e.*) = ^($-$-)^-^C0s(2.^+^). 

(4.1) 

Here is the created by the external field flux thread- 
ing through the loop, ^ = LI and the phase shift S 
is equal to and tt for 0- and 7r-junction, respectively. 
The first term in Eq. (4.1) corresponds to the screening 
whereas the second term is standard for an Josephson 
junction. The important question is if the spontaneous 
magnetic moment (fiux or supercurrent) occurs in zero 
external field. The answer to this question may be ob- 
tained by finding minimum of energy (4.1). 
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FIG. 19. Flux dependence of the dimensionless free energy 
F of the one-loop model. For the 0-junction with L = 20 
(upper panel) and 7r-junction with L — 0.5 < 1 (dotted line 
in lower panel) the function Fhas minimum at $/$o = 0. 
For a TT-junction with L = 20 > 1 (solid line in lower panel) 
the minimum is located at $/3>o and the paramagnetic 
response becomes possible. 
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In the absence of the external field the dimensionless 
energy F($/$o) reads as 



F($/$o) = 



27r2 



cos 



2L 



$0 



(4.2) 



Here the dimensionless self-inductance L is equal to 



(4.3) 



It is easy to show that at $/$o = the energy has ex- 
tremum and its second derivative with respect to $/$o 
has the simple form 



^"'(0) = 2(1 ±L) 



(4.4) 



where the sign phis and minus corresponds to the 0- 
and TT-junction, respectively. For an unfrustrated loop 
-F"(0) > and the state without the supercurrent should 
be stable and the PME is not, as expected, possible. 
In the case of the TT-junction the energy has minimum 
$/$o = (F"(0) > 0) for Z < 1 and maximum{F" (0) < 
0) for L > 1. We come to a very interesting result: the 
energy of a frustrated loop has minimum at <l> ^ pro- 
vided Z > 1. It is demonstrated in Fig. 19. Thus, above 
the critical value Lc = 1 the spontaneous magnetic mo- 
ment Msp appears in the single TT-junction loop leading 
to the PME. 

In the L — > oo limit Mgp may be determined ex- 
actly. In this case the second term in (4.1) dominates 
and the energy has a ladder of minima at $ = n$o and 
<I> = (n -|- i) $0 is integer number) for 0-junction and 
TT-junction, respectively. The spontaneous magnetization 
corresponding to the global minimum with n = is equal 
to Msp = ^o/8ttS, where S is the loop area. 

Since the energy minima of the single TT-junction are 
located near 

$ = (n -t- i) $0 the PME IS, m some sense, 
the manifestation of the common half-integer flux quan- 
tum effect [15]. The screening plays a key role in the 
observation of the paramagnetic response. Taking into 
account the temperature dependence of critical current 
Ic and using the one-loop model (4.1), Sigrist and Rice 
were able to reproduce the paramagnetic behavior in the 
FC regime [9,10]. In accord with experiments the para- 
magnetic signal is suppressed as the external field is in- 
creased. 



B. Multi-loop model 

Although the single-loop model considered in the pre- 
vious chapter predicts the appearance of the PME, it fails 
to capture many cooperative phenomena such as the ag- 
ing, the compensation effect, a possibility of the existence 



of a new chiral glass phase etc in granular materials. This 
shortcoming may be amended in the multi-loop model 
[12,145] where the interaction between loops is taken into 
account. 

B.l. Hamiltonian 

Suppose that weak links connecting the neighboring 
grains are distributed sufficiently dense, so that the sys- 
tem can be viewed as an infinite network of Josephson 
junctions which are not decomposed into finite clusters. 
We model such ceramic superconductors by a hypercubic 
(two- or three-dimensional) lattice model of a Josephson 
junction arrays with finite self-inductance. Neglecting 
charging effects of the grain the Hamiltonian is given by 
[12] 



cos I 



<ij> 



2^ <V> *° 



where 4>i is the phase of the condensate of the grain at 
the i-th site of a simple hypercubic lattice, A is the fluc- 
tuating gauge potential at each link of the lattice, de- 
notes the Josephson coupling between the i-th and j-th 
grains, and L is the self-inductance of a loop. The effect 
of screening currents inside grains is not considered ex- 
plicitly, since for large length scales they simply lead to a 
Hamiltonian Ti with an effective self-inductance L [146]. 
The first sum is taken over all nearest-neighbor pairs and 
the second sum is taken over all elementary plaquettes on 
the lattice. Fluctuating variables to be summed over are 
the phase variables, (f)i, at each site and the gauge vari- 
ables, Aij, at each link. $p is the total magnetic flux 
threading through the p-th plaquette, whereas $p^* is 
the flux duo to an external magnetic field applied along 
the 2-direction, 



ext _ j HS if p is on the < xy > plane 
I otherwise , 



(4.6) 



where S denotes the area of an elementary plaquette. 
The Hamiltonian (4.5) is an extension of the single-loop 
energy (4.1) to the interacting loops case. 

In the s-wave ceramics, sign of the Josephson coupling 
is always positive (Jy > 0), while in ci-wave ceramics 
it could be either positive (0-junction) or negative (tt- 
j unction) depending on the relative direction of the junc- 
tion and the crystal grains on both sides. In d-wave case, 
the sign of Jij is expected to appear randomly since the 
spatial orientation of each crystal grain would be random. 
It was also suggested [147] that the Josephson junction 
between d-wave superconductors under certain circum- 
stances (e.g., near the interface) might have an energy 
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minimum at some fractional value, neither at or tt, as a 
result of a spontaneous time-reversal symmetry breaking 
at the junction. Such a possibility is not discussed in this 
review. 

The model with uniform ferromagnetic couplings 
{Jij = J) studied by Dasgupta and Halperin [148] 
exhibits a standard normal-superconductor transition. 
Here we deal with two types of bond distributions. For 
d-wave ceramics Jij is assumed to take the value - J (tt- 
junction) with probability c and +J (0-junction) with 
probability (1-c) [12]. In the s-wave model Jij is always 
positive ("ferromagnetic" interaction) but distributed 
uniformly between and 2tt [145]. 

In the case of d-wave superconductors the frustration 
due to competition between "ferromagnetic" and "an- 
tiferromagnctic" interactions [144] causes the nonzero 
spontaneous magnetic moment provided c exceeds the 
percolation threshold. The global PME becomes, there- 
fore, possible. In the absence of an external field, the 
global minimum of the energy of the s-wave system corre- 
sponds to the configuration with zero supercurrents and 
the system displays a standard diamagnetism. 

B.2. Critical self-inductance 

In the single-loop model with one 7r-junction the di- 
mensionless critical inductance above which the system 
shows the PME is equal to L = 1. The question we 
as now is does the interaction between frustrated loops 
change this critical value. For the multi-loop model (4.5) 
L is also given by Eq. (4.3) with the critical current [145] 



27rJ 
$0 



(4.7) 
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FIG. 20. The temperature dependence of the FC suscepti- 
bility (a) and of the ZFC susceptibility (b) of d-wave ceramics 
for several values of the external field. The lattice size is equal 
to Z = 8, L = 1 and the concentration of n junctions c = 0.5. 
After Kawamura and Li [145]. 

B.3. Simulation results 



One can show (see Appendix A) that, contrary to the 
one-loop model, the critical inductance above which the 
PME is observable is equal to 



Total magnetization along the z axis normalized per 
plaquette m is given by [12] 



0. 



(4.6 



In the limit of small inductances the spontaneous flux 
of the multi-loop model is equal to [145] 



$0 



— RJ ± 



for L < 1 . 



(4.9) 



m = 



(4.10) 



pE<xy> 



where the sum is taken over all Np plaquettcs on the 
< xy > plane. The linear susceptibility x (dimensionless 
in CGS units) may be calculated via fluctuation of the 
magnetization as follows 



The last expression shows that for small inductance, 
spontaneous flux remains nonzero but becomes small. In 
real ceramic samples one has an additional diamagnetic 
contribution from intragranular supercurrents neglected 
in the present model, a small paramagnetic contribu- 
tion from the intergranular supercurrents may easily be 
masked. Therefore, in practice, moderately large induc- 
tance would be necessary to observe the paramagnetic 
behavior. 



X = -— f[<m2>-<m>2]j-— , (4.11) 
ksTL 47r 



where < ... > and [...] j represent a thermal average and 
a configurational average over the bond distribution, re- 
spectively. For conventional superconductors the first 
term in (4.11) vanishes at low temperatures and x takes 
on a standard value 
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FIG. 21. The temperature dependence of the ZFC suscepti- 
bility (a) and of the ZFC susceptibility (b) of s-wave ceramics 
for several values of the external field. The lattice size is equal 
to Z = 8 and L = 1. After Kawamura and Li [145]. 

The temperature and field dependence of the FC and 
ZFC susceptibility was studied by Langevin dynamics 
[12] and by Monte Carlo simulations [145]. Both these 
methods give essentially the same results and we present 
those obtained by the latter approach. 

As can be seen from Fig. 20, a paramagnetic behavior 
is shows up in FC regime in the field range h <1 (roughly 
corresponding to ~ $o/4), where the dimensionless 
field h is defined as [145] 



h = 



2ttSH 

$0 



(4.12) 



By contrast, xzfc remains diamagnetic at low temper- 
atures for any value of the dc fields studied. 

From Fig. 20 it follows that the remanent magnetiza- 
tion nirem ~ ™FC — mzFC is positive, consistent with 
experiments on high-Tg ceramics [149,150]. This obser- 
vation indicates the occurrence of flux trapping in the 
d-wave model. 

The Monte Carlo simulations [145] showed that the 
PME is clearly visible for a very small inductance L = 
0.1. This result compatible with the conclusion in the 
previous section that Lc = for the multi-loop model. 



Fig. 21 shows the temperature and field dependence of 
Xfc and xzfc of the s-wave system. We set the system 
size I = 8 and L = 1. As expected, the response is al- 
ways diamagnetic. Furthermore, there is no appreciable 
difference between xfc and xzfc for h < 0.2. Within 
the statistical errors rrirm turns out to be zero confirming 
that the flux does not get trapped inside a sample 



C. Other models 

Using a Josephson junction model of special geometry 
[151] Auletta et al. [152] have shown that the PME is 
possible even without 7r-j unctions. The similar conclu- 
sion was also made by other groups [153,46,47,154,155] 
for standard systems of 0-junctions (the effect of mutual 
inductances has been considered in [153,155]). The para- 
magnetic behavior can arise either in a d-wave sample 
or in a conventional sample of confined geometry. From 
this point of view, the PME observed in simulations of 
the 0-junction networks [152,153,46,47,154,155] may be 
due to the finite size effect or it is merely a metastable 
effect. 

Recently, Khalil [20] assumed that near randomly 
distributed defects the localized superconducting states 
with different orbital moments could exist. The first- 
order transitions between these states below the critical 
temperature are accompanied by jumps of the magnetic 
moments. As a result, multiple-quanta vortices and spon- 
taneous moments randomly distributed in the material 
matrix will be formed. This provides the explanation of 
the PME. 

Authors of Ref. [156] proposed that a certain amount of 

local moments exist that are becoming partially aligned 
by an external field. Based on experimental data, these 
moments are estimated to be too small compared to the 
half flux quanta. They are can not be, therefore, due to 
TT- junctions but their nature remains unspecified. 

The another interpretation [157] of the PME is based 
on the so called impurity mechanism [158] for high- Tc su- 
perconductivity. The validity of this heuristic argument 
remains ambiguous. 

Finally, the PME was proposed to be related to the 
vortex pinning [2] at the Kosterlitz-Thouless transition. 
Since this transition is possible in two dimensions [159] it 
is not clear if the scenario of Svelindh et al [2] is applied 
to three-dimensional systems. 



V. FLUX COMPRESSION AND 
PARAMAGNETIC MEISSNER EFFECT 

In this chapter we discuss the fiux compression mech- 
anism for the PME which is probably applied to conven- 
tional superconductors with s-wave type of pairing sym- 
metry. There are two approaches to study the flux trap- 
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ping in the confined geometry. One of them is based on 
the Bean critical model and was developed by Koshelev 
and Larkin [19]. The second approach [21] relied on the 
numerical solution of the Ginzburg-Landay equations, 
predicts the paramagnetic response due to spontaneous 
magnetic moments in giant vortex states which occur be- 
low the third critical field. 



A. Bean model 

The PME observed in conventional superconductors 
[18] and in twined single crystals of YBa2Cu307_5 sug- 
gests the possibility of alternative explanations without 
exotic d-wave pairing. One of such possibilities is based 
on the flux compression mechanism [19] using the Bean 
model [160] for the critical state . A possible origin of 
the existence of the spontaneous paramagnetic moment 
is flux trap caused by an inhomogeneous superconduct- 
ing transition. If the edges of the sample becomes su- 
perconducting first due to, e.g., inhomogeneity of mate- 
rial (surface layers have a higher value of Tc compared 
to the bulk) or inhomogeneous cooling, then expulsion 
of vortices from this area opens a flux-free region near 
the edges. This gives rise to the flux capture inside the 
sample. Further cooling leads to a broading of the flux- 
free region and further compression of the flux due to, 
e.g., the vortex Nernst effect [161]. As the whole sam- 
ple becomes superconducting the critical Bean state de- 
velops in the flux region. For a thin superconducting 
strip geometry (see Fig. 22) one can imagine this sce- 
nario as follows. The outer current mimics the diamag- 
nctic shielding current, while the inner current models 
the paramagnetic pinning current. If the fltix is trapped 
inside sample, then the flux created by the outer current 
should be completely compensated by the flux generated 
by the inner current in the flux free region. Since the field 
generated by currents flowing in the inner region changes 
sign, these currents give a smaller contribution to the to- 
tal flux than do the currents flowing in the outer region 
[19]. This should be compensated by a larger magnitude 
of the inner paramagnetic currents, which gives rise to 
the paramagnetic moment. 

We now give a qualitative description of the existence 
of the PME due to flux compression in a simple case of 
a thin superconducting slab of width 2w and thickness 
d situated in a transverse magnetic field [19] as shown 
in Fig. 22. The Bean state is assumed to be complete 
with the critical current jc fiowing in the strip of width 
2b {\x\ < b < w) surrounded by flux free region of width 



H{z) = - 



Hj 



1/2 



In 



w ■ 



w — b (b < \x\ < w) (the incomplete Bean state was also 
discussed by Koshelev and Larkin [19]). 



Comi:le*ti? RS'^in state 




FIG. 22. Complete Bean critical state in a thin supercon- 
ducting slab after field cooling. Reprinted from Koshelev and 
Larkin [19]. 

Then the tangential (H^) and longitudinal (Hy) field 
components satisfy the following boundary conditions 
[19] 



Hy 



= Hjsign{x), < |a;] < b, 
= 0, \x\ > b, 

= 0, b <\x\ < w, 



00 



(5.1) 



where 2Hj = (Att/cJJc, Jc = djc- The first equation 
in (5.1) from the Gaussian theorem and the continuous 
boundary conditions for the tangential component of the 
field near the slab surface [162,163]. The second and third 
equations are valid for the flux-free region. The field dis- 
tribution may be obtained from the Maxwell equations 



divB = , rotB = 



(5.2) 



subjected to boundary conditions (5.1. An elegant way 
to solve these equations is based on theory of complex 
variable functions [19]. Introducing the complex function 
H(z) — Hx + iHy of complex variable z = x + iy, then 
from Maxwell equations (5.2) it follows that Hx and Hy 
satisfy the Cauchy conditions for an analytic function. 
The problem of finding H{z) is, therefore, reduced to the 
finding of an analytic function satisfying conditions (5.1). 
One can show that [19] 



— In 



{i-zybY^^ + {i-zyw^)^/^ 

(1-^762)1/2 -(1-^2/^2)1/2 



H^ 



1/2 



(5.3) 



Using this function we obtain the total magnetic fiux 



19 



.6 

! / Hydx 
Jo 



-^^^{[E(m) - (1 - m)ii:(TO)]/ni^ — ^ 
TT 1 — vm 
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2y/mK{m)} 



(5.4) 



where E{m) and K{m) are complete elliptic integrals 
[164]. The condition $ = 2wH^f, with / being the 
fraction of trapped flux, gives the relation between the 
critical field Hj and the width b of flux region. 

The magnetic moment can be determined from the 
large-distance asymptotics of the field as 



Hiz) 



which gives 
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(5.5) 
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(5.6) 



Parameter m has to be expressed through the critical 
current Jc = cHj /2TT and the fraction of captured flux 
/ = <^/2wH^ using Eq. (5.4). In the Hmit to << 1 
(strong compression) the magnetic moment 



M = 



Hjw^/m (A 1 

~ 7 



+ 0(to) 



(5.7) 



which is positive for / > 7r/4 « 0.79. Remarkably, this 
value almost coincides with the critical value / « 0.8 
(above which a paramagnetic moment is possible) ob- 
tained numerically [19]. Fig. (23) shows the dependence 
of M on the relative width of flux region h/w for /=1 
obtained from Eq. (5.4) and (5.6). 



< a I * 



1 



FIG. 23. Dependence of magnetic moment in the com- 
pressed flux state (/ = 1) on the relative width of the flux 

region. For comparison similar dependence for an ideal su- 
perconductor {jc = 0) is shown. After Koshelev and Larkin 
[19]. 

In the weak compression regime w — b « w one can 
show that the magnetization is always negative [19] . 



B. Ginzburg-Landau equation approach 

In this section we consider the flux compression mech- 
anism by the another approach [21] which is based on 
the self-consistent numerical solution of the Ginzburg- 
Landau equations. One can show show that due to con- 
flned geometry the giant vortex state is formed in the FC 
regime at the third critical field HcsiT). The PME ap- 
pears as a result of the flux trap inside the giant vortex 
state with a fixed orbital quantum number L. 

The basic postulate of the Ginzburg-Landau theory is 
that if the Cooper wave function ^ is small and varies 
slowly in space, the free energy can be expanded in a 
series of the form [165] 



2to* 



i c 




(5.8) 



where a{T) cx (T—Tco), Tco is the mean field transition 
temperature. Minimizing the Ginzburg-Landau free en- 
ergy with respect to the fields ^'(r) and vector potential 
A{r) one can get two well-known equations which will be 



used to study the PME. 

Following Ref. [21] (see also Ref. [166]) we consider 
a infinitely long cylinder of radius ro in a magnetic field 
parallel to its axis. First we will find Hcz{T) above which 
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the miclcation of the superconducting phase takes place 
and then study the PME. Determination of Hc^{T) re- 
quires to solve the linearized Ginzburg-Landau equation 
for the superconducting order parameter [7] 



1 
2m 



(5.9) 



which is identical to the Schrodingcr equation for a par- 
ticle of charge e* = 2e in a uniform magnetic field given 
by the vector potential ^ = V x ^4. Parameter a in 
Eq. (5.9) plays the role of energy E in the corresponding 
Schrodinger equation and 



E 



T, 



cO 



(5.10) 



where ^(T) is the bulk coherence length and Tgo is the 
critical temperature at zero field. 

For finite samples the order parameter 5* should satisfy 
the boundary condition for a superconductor-insulator 
interface [167] 



-inv ■ 



^-L = 



(5.11) 



which assures that no current passes through the surface. 

The common strategy to find HcsiT) is that [7] one has 
to solve the linearized Ginzburg-Landau equation sub- 
jected to the boundary condition (5.11). Then Hc^iT) 
is subtracted from the lowest energy E{H) which gives 
the highest T in Eq. (5.10) coinciding with the nucle- 
ation phase boundary. In the case of large bulk samples 
when the boundary condition (5.11) becomes irrelevant, 
the lowest Landau level E = TiLoj^ (the cyclotron fre- 
quency io = e*H/mc) gives the well known upper critical 
field Hc2 



he T, 



cO 



$0 



c2 



277^2 (T) ■ 



(5.12) 



For samples of confined geometry the boundary con- 
dition (5.11) must be taken into account. The crucial 
question is how small the sample should be to observe 
the finite size effect. To answer this question one may 
refer to analogy between the Ginburg-Landau equation 
(5.9) and the Schrodinger equation for normal electron 
considered by Dingle [168] for the analysis of quantum 
oscillations in small metallic samples. The finite size ef- 
fect is expected to become visible [168] provided ro < r^, 
where ri, is the Larmor radius or 



H X ro < 5{Gcm) . 



(5.13) 



Clearly, the critical size depends on the applied mag- 
netic field. If the size is equal to, say, ro = 5^m then 
the boundary condition (5.11) must be taken account for 
H < 10* G. 



In order to solve the linearized Ginzburg-Landau equa- 
tion (5.9) for a long cylinder we choose the cylindrical 
coordinates (r, (p, z) and the following gauge 

A = (0,i/r/2,0) . (5.14) 
Then the solution reads as [169] 
^LiRA) = e=^^^'^i?^exp(-i?V2)M(-7V,L + l,i?2) 

(5.15) 

and the energy E± of motion in the plane perpendicu- 
lar to H is given by the orbital quantum number L and 
parameter N: 

= ^^{2N±L + L + 1) . (5.16) 
2mc 

Here M{-N, L + l,R'^)is the Kummer function [169], the 
dimensionless radius R = ^/7r, 7 = e* H / {2hc). From the 
condition $(r, ^ + 27r) = $(r, (j)) it follows that L should 
be integer number. It is important to stress that N is not 
necessarily integer number and it must be found from the 
boundary condition (5.11) which with the gauge choice 
(5.14) is reduced to 



d\m)\, 

dR 



I R=Rq 



= 0. 



(5.17) 



So far as we are looking for the lowest possible energy 

state for determination of Hc?,{T), one should take the 
minus sign in the argument of the exponent exp(—iL(j)) 
in the solution (5.15). Then the energy levels become 



Ej_ = huj{N+ 1/2) 



(5.18) 



where u) is the cyclotron frequency. This expression co- 
incides with the well-known Landau levels but now N is 
any real number including negative one. 

From Eqs. (5.15) and (5.17) we obtain the equation 
for determination of the parameter A'^ [21] 



{L^ - Rl)M{-N,L + l,Rl) 



^M(-iV. 



l,L + 2,Rl) = 0. 



(5.19) 



Remarkably, the numerical solution of the last equation 
gives negative values for N{L, R.o) which immediately 
lead to the energy E± in (5.18) lower than the bulk value 
huj/2. Therefore, as seen from Fig. 24, in the presence 
of the finite surface the superconductivity can occur well 
above Hc2{T) line obtained for A'' = 0. The cusplike 
phase boundary H^s occurs due to switching between dif- 
ferent orbital momenta L. Such a phase boundary was 
observed experimentally [170,171] for the superconduct- 
ing disks. The linear component of the cusplike Hcs line 
is 1.965i?c2, which is in good agreement with the calcu- 
lations of Hc3 in the L — > 00 limit [167,172]. 
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In order to study the PME one has to answer the 
following question: how to mimic the FC and ZFC ex- 
periments by the Ginzburg-Landau method? In the FC 
mode when the external field is applied to a sample the 
states with orbital quantum numbers L > correspond- 
ing to the rotation of the superconducting condensate 
due to the action of the Lorcntz force arc expected to 
be formed. Then by cooling down we are crossing the 
Hc3{T) boundary (see Fig. 24) at a given L > and 
reaching a temperature T < Tcz{H) at which measure- 
ments are carried out. It should be noted that [173] L 
can be retained over a large temperature interval, namely, 
from T{Hc7i) to T{Hc) (near or at T{Hc) the quantum 
number L changes abruptly, thereby expelling a large 
amount of flux from the sample, as the specimen makes 
a transition to the Meissncr state). From the physical 
point of view, the conservation of L is related to pin- 
ning of the giant vortex state by the sample boundary. 
Thus, the FC experiments may be analysed by solving 
the Ginzburg-Landau equations below Hc^iT) for fixed 
L > Q. The same is valid for the ZFC regime but with 
L = Q. 




''L=3 




FIG. 24. The third critical field iJcS and the bulk upper 
critical field Hc2 (dashed line) vs normalized temperature 
t = T/Tc. The cusplike Hc3 is formed due to the jump of 
the orbital quantum number L. The three-dimensional plot 
of 1*^1 is shown for several values of L. After Moshchalkov et 
al. [21]. 



Below Hc3{T), one needs to solve the full Ginzgurg- 
Landau equations because the superconducting conden- 
sate induces the screening due to supercurrents (propor- 



tional to which should be taken into account. Then 
two coupled dimensionless equations are [174] 

(^-^ * = *(1- 1*1^) , 
V X V X A = ^ 




(5.20) 



where and r are in units of ^infty,V^HcX{0) and 
A(0), respectively, 'i'oc, He, A(0) and k arc the wave bmc- 
tion in bulk, the thermodynamic critical field, the pen- 
etration depth at zero field and zero temperature, and 
the Ginzburg-Landau parameter, respectively. For the 
cylindrical symmetry we choose 



*(r,(?i) = e'^'^F(r), A= (0, 



Hr 



2r 



,0). (5.21) 



Comparing Eqs. (5.21) and ((5.14) one can see that the 
vector potential in the full Ginzburg-Landau equations 
contains the induced field neglected in the linear approx- 
imation. Denoting the local magnetic field by b and using 
6 = V X ^ we have 



H 



1 dip 
2r dr 



(5.22) 



Substitution of Eqs. (5.21) and (5.22) into Eq. (5.20) 
gives two coupled equations for determination of F and 
if [21] 



dr^ 



IdF _^ 
r dr 



1 1 Kip L 
-KHr+-— 

2 2 r r 



fip 
dr^ 



1 d(p 
r dr 



+ {Hr'+ip 



F - kF{1 - F'^) 
2L\ 



(5.23) 



The corresponding boundary conditions are [175] 



dF _ ^ dip 
dr ' dr 



= 



(5.24) 



at r = ro and 



ip = 0, F = for L 7^ , 
dF 

^ = for i = 
dr 



(5.25) 



at r = 0. 

The magnetization per unit volume 47rM can be de- 
fined via the flux through the cylinder and has the fol- 
lowing form 



AttM 



pro 

I (Kr) 
Jo 



H)rdr = ip{ro) ■ (5.26) 
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nctic response corresponding to a transition from large L 
values trapped at ffcs to the state with L = 0. 

The numerical solution of Eqs. (5.23) - (5.25) for the 
ZFC regime with L = gives, in agreement with experi- 
ments, the diamagnetic response at any temperature [21]. 
Thus, the PME and its field dependence can be obtained 
from the self-consistent solution of the Ginzburg-Landau 
equations assuming that orbital quantum number L is 
conserved. This very simple and natural approach does 
not involve any further assumptions related to the exis- 
tence of a TT junction or d- wave superconductivity. The 
drawback of the Ginzburg-Landau equation approach is 
that it docs not allow for studying the PME far from the 
transition temperature to the superconducting phase. 

The PME in a mesoscopic disk has been also studied 
[176,55] by the Ginzburg-Landau method . The results 
are qualitatively the same as for the case of the cylindri- 
cal symmetry. 

To summarize this section, according to the Ginzburg- 
Landau theory, the paramagnetism in low-Tg mesoscopic 
superconductors may be caused by moments in giant vor- 
tex states. The vortices formed inside a sample in the 
FC regime get pinned to the boundary which is a source 
of inhomogeneity. The directions of the external and of 
the vortex magnetic fields are the same but the direc- 
tions of currents screening these fields from the bulk are 
opposite. Since the screening currents contribute to the 
magnetic moment they give two contributions of opposite 
sign (M = M^ia + Mpara)- The distributions of the fields 
and currents inside the sample vary with H and the mea- 
sured magnetization, therefore, may be either positive or 
negative depending on H. 

VI. CHIRAL GLASS 



FIG. 25. Temperature dependence of Mfc of a supercon- 
ductor with the Ginzburg-Landau parameter k = 10 for dif- 
ferent values of L in applied magnetic field H = Q.l (a) and 
H = 0.01 (b). The radius of the cylinder ro = V3A(0). After 
Moshchalkov et al. [21]. 

As in the linear approximation the orbital quantum 
number L remains a good quantum number due to the 
cylindrical symmetry. For a fixed L one can solve Eqs. 
(5.23) - (5.25) numerically to obtain the magnetization 
in the FC mode. The typical results are shown in Fig. 
25. In accord with experiments the PME occurs for the 
fixed values of L and it gets enhanced as the external 
field is lowered. The paramagnetic signal is, however, 
sensitive to H and L. For L = 1 and 2 the increase in 
the field from H = 0.001 to H = 0.001 changes the sign 
of the magnetization from unusual paramagnetic to con- 
ventional diamagnetic. The reported sensitivity to the 
surface treatment [16,17,42] may be caused by violation 
of the L conservation and recovery of a normal diamag- 



The existence of spontaneoiis superciirrents in ceramic 
superconductors may lead to, as assumed by Kusmartsev 
[8] , a so called " orbital glass" . A natural question arises 
is what is the nature of the orbital glass? The picture by 
Sigrist and Rice [9] is essentially a single-loop picture in 
which the interactions between loops are irrelevant. In 
their scenario, when the temperature is lowered across 
the superconducting transition of each grain T/*" , the FC 
susceptibility changes sign from negative to positive at a 
certain temperature Tq, sightly below T^''. Such a change 
of sign of xfc , which is often regarded as a measure of the 
orbital glass transition point, is a crossover not related to 
any intergranular (interloop) cooperative phenomena. In 
this picture, the PME arises as a property of an ensemble 
of noninteracting loops if there occurs an intragranular 
superconducting transition. By contrast, the possible co- 
operative character of the orbital glass was pointed out 
by Dominguez et al [12] and by Khomskii [6], although 
these authors did not detail the nature of the cooperative 
phenomena. 
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Naturally, one may ask here whether could be a ther- 
modynamically stable orbital glass state characterized by 
a spontaneous breaking of certain symmetry over an en- 
tire granular system. An issue to be addressed is whether 
there could be some sort of thermodynamic intergranu- 
lar phase transition accompanied with a divergent length 
scale over grains, and if it is, what is the order parameter 
of such phase transition. It is reasonable to assume that 
the orbital glass is, similar to spin glass, a state of frozen 
spontaneous moments or fluxes due to frustration caused 
by random distribution of tt junctions. 

An attempt [177] to detect such a orbital (or flux) glass 
by Monte Carlo simulations with the use of multi-loop 
model (4.5) gave an ambiguous conclusion. Instead, the 
simulations [27,177] showed clearly that a novel thermo- 
dynamic phase may occur in zero external field in certain 
ceramic high-Tg superconductors. This phase is char- 
acterized by a spontaneously broken time-reversal sym- 
metry with keeping U{1) gauge symmetry, and is called 
"chiral glass phase". The order parameter is then a "chi- 
rality" , quenched-in half a vortex, which represents the 
direction of the local loop-supercurrent circulating over 
grains. As one can see below, since chirality and flux 
have the same symmetry we believe that the orbital glass 
should also exist in three dimensions. The fact that it was 
not identified in Monte Carlo simulations [177] may be 
related to small system sizes. 

In this chapter we introduce the chirality concept in 
discuss the chiral glass phase in detail. Experimental 
search for this new phase will be presented. 



A. Chirality concept 



respectively. One may define the chirality n via a vec- 
tor product of two neighboring spins averaged over three 
spin pairs 

<ij> <ij> 

(6.2) 

Clearly, k = 1 for the right-handed (clockwise) configu- 
ration in left panel of Fig. 26 and k = —1 for the left- 
handed (counterclockwise) configuration in right panel. 
In the case of XY spins considered here, the chirality k 
is actually a pseudo scalar. It remains invariant under 
global S0{2) proper spin rotations while it changes sign 
under global Z2 spin reflections. In order to transform 
the chiral state with k = 1 to the state with k = — 1 one 
needs to make a global spin reflection. 

It should be noted that if the interaction is ferro- 
magnetic then the system becomes geometrically unfrus- 
trated (all spin in the ground state are parallel) and 
K = 0. So the chirality concept makes sense only for 
frustrated systems. The chirality defined by Eq. (6.2) 
is a vector quantity, while in the literature one discusses 
also a scalar chirality [178] for Heisenberg spins. 




Frustration in vector spin systems often give rise to 
noncollinear and noncoplanar spin orderings which may 
be characterized by the so call chiral degrees of freedom. 
While the chirality concept has long been a familiar con- 
cept in molecular chemistry, it was introduced into the 
field of magnetism first by Villain [28] . 

The simplest way to introduce the chirality concept 
[28] is to consider the XY model defined by the following 
Hamiltonian 

H = JijSi.Sj = Jij cos{4>i - 4)j) , (6.1) 

ij ij 

where the unit spin vector Si = (cos (/>j, sin (/>j). Model 
(6.1) is the simplest version of (4.5) without the external 
field and screening. We consider the triangular lattice 
and assume that the spin-spin interaction is uniform and 
antiferromagnetic (Jjj = J < 0), then the ground state 
of three spins located at corners of a triangle, as shown in 
Fig. 26 is two-fold degenerate. Angles between neighbor- 
ing spins are 120° and -120° for the left and right panel. 



FIG. 26. The ground state spin configuration of three 
XY spins antiferromagnetically coupled on a triangle lat- 
tice. Frustration leads to the canted ordered state. The 
ground state is twofold degenerate according to whether the 
non-coUinear spin structure is right- or left-handed, each of 
which is characterized by the opposite chirality as shown in 
left and right panel. 

The chirality concept plays an important role in 
physics of classical [178,179] and quantum [180] spin sys- 
tems. Especially since the discovery of high- Tc cuprates 
the scalar chirality has been a key theoretical concept 
in physics of strongly correlated electrons [181-183]. In 
order for the spin chirality to be ordered both the time- 
reversal and parity symmetry must be broken. Break- 
ing of these symmetries in two dimensions brings about 
many intriguing physics such as parity anomaly [184,185] , 
anyon superconductivity [186] , and quantized Hall effect 
in zero external magnetic field [187]. In this review we 
try to attract the attention of readers to a possibility to 
observe the so called chiral glass phase in ceramic mate- 
rials. 
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B. Simulation search for the chiral glass phase 

In order to study the nature of ordering of the chiral 
glass phase in ceramic superconductors wc use the multi- 
loop Hamiltonian (4.5) with the zero external field. The 
Josephson coupling Jij is assumed to be an independent 
variable taking the values J or — J with equal probabil- 
ity (bimodal distribution). While we present results for 
this particular distribution, one could expect from ex- 
perience in spin glass studies that the results would be 
rather insensitive to the details of the distribution. 

Contrary to the gauge-glass model [188-190], the 
Hamiltonian (4.5) defined at i7 = 0, keeps the Z2 time- 
reversal symmetry in addition to the U{1) gauge symme- 
try. Frustration arises from the random distribution of 
both positive and negative Josephson couplings. This 
should be contrasted to the vortex-glass (gauge-glass) 
problem, where the associated Hamiltonian does not pos- 
sess the time-reversal symmetry due to external magnetic 
fields, while the frustration arises from the magnetic field 
but not from the Jij . 

Extending the definition (6.2) to the case of the multi- 
loop model the local chirality at each plaquettes is given 
by the gauge- invariant quantity [145] 

p 

Kp = 2"^/^ ^ sign( Jy ) sin((/)j - - Aij), (6.3) 

<ij> 

where the sum runs over a directed contour along the 
sides of the plaquette p. Physically, the chirality, Kp, 
is a half (tt) vortex, being proportional to the loop- 
supercurrent circulating round a plaquette p. If the pla- 
quette p is frustrated, the local chirality Kp tends to 
take a finite value, its sign representing either clockwise 
or counterclockwise orientation of circulating supercur- 
rent. If, on the other hand, the plaquette is unfrus- 
trated, as mentioned before, it tends to take a value 
around zero. Note that the chirality is a pseudoscalar 
in the sense that it is invariant under global L''(l) gauge 
transformation, ^ + A(j), , but changes 

its sign under global Z2 time-reversal transformation, 
^ij ~^ Due to this symmetry prop- 

erty, chirality can be regarded as an order parameter of 
the chiral order. 

In analogy with the spin glass theory, the chiral glass 
phase may be characterized by the Edwards-Anderson 
order parameter [191] Qcg- 



4^ 



[<«p>']j. (6.4) 

In the chiral glass phase c^q ^ 0, while in the chirality- 
disordered ("paramagnetic") phase q^Q — 0. In the chi- 
ral glass state chiralities are orientated randomly in space 
but frozen in time. 

The magnetization and linear susceptibility are given 
by Eqs. (4.10) and (4.11) but with = 0. The non- 
linear susceptibility, X2 is defined as follows [145,177] 



X2 



1 (Pi 



1 /ttJTV, 



6di73 



4 < TO >< TO^ > -I- 



-3 < TO^ >2 

12 < X m -6 < m >^lj 



(6.5) 



Note that X2j being proportional to the minus of the 
third-harmonic component of the ac susceptibility, is 
sometimes denoted as X3 in the literature. 

Following the standard technique from the spin glass 
theory [144] we introduce the overlap q between the chiral 
variables in the two independent replicas [27,177] 



(6.6) 



where upper indices 1 and 2 denote two replicas. In terms 
of this chiral overlap the Binder ratio of the chirality is 
calculated by 



9CG 



3- 



[< >]./ 
[< q^ >H 



(6.7) 



Here gcG is normalized so that in the thermodynamic 
limit it tends to zero above the chiral-glass transition 
temperature. Tog, and tends to unity below Tcg pro- 
vided the ground state is non-degenerate. At the chiral- 
glass transition point, curves of gcG against T for dif- 
ferent system sizes / should intersect asymptotically. It 
should be noted that one can use not only the standard 
Binder function g but also so called A and G function 
[192-194] to study the nature of ordering in frustrated 
systems. 

The chiral-glass susceptibility, which is expected to di- 
verge at the chiral-glass transition point, is given by 



XCG = Np[< q^ >]j . 



(6.8) 



In the framework of the one-parametric scaling theory 
the correlation length ^ diverges as one approaches to the 
chiral glass transition temperature by a power law 



\T-Tcg\ 



(6.9) 



The divergence of the nonlinear and the chiral glass sus- 
ceptibilities near Tcg is characterized by critical expo- 
nent 72 and 7cG 



X2 ~ \T-Tcg\-^' , 
XCG ~ |T-Tcg|"^^° 



(6.10) 



For the finite size scaling analysis we introduce the 
scaling functions [144,27,177] 



gcG = gcGil^'"'"' \T-TcG 
for the Binder parameter. 



(6.11) 
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X2 = L^^/''^°x2(/'/"^° \ T-TcG I), (6.12) 
for the nonlinear susceptibility and 

XCG = L^-^^^xcGil'/"'"' \ T-TcG I) (6.13) 

for the chiral glass susceptibility. The critical exponents 
can be determined from the condition that the scaling 
functions, plotted versus argument I'^^'^cc | y _ j"^^ |^ 
should not depend on system size I [195]. In the other 
words, they must be collapsed onto a single curve. Expo- 
nent riGG determining the decay of the chiral correlation 
function at T = Tea ([< Kp(0)Kp(r) >],/ - r-^+^-'"^<^, 
where D is spatial dimensionality) and exponent jcG are 
related by a simple scaling law 



ij.4 
L.J 

c.a 

L.l 

r ci 

C.i 




7CG 



(2 - 7]cg)vcg 



(6.14) 
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In the vicinity of the transition point the Edwards- 
Anderson order parameter behaves as 

qE^ r^lT-TcGf'"' . (6.15) 

Exponent pcG is expressed via the other exponents by 

2/?CG = DvcG - ICG ■ (6.16) 

It should be noted that the chiral glass phase in the 
multi-loop model (4.5) without screening has been stud- 
ied [196,197] by Monte Carlo simulations. It was shown 
that thermal fluctuations destroy the three-dimensional 
spin glass ordering but leaving the chiral glass ordering 
to be stable at nonzero temperature. Although chirali- 
ties and Ising spins have the same symmetry, the chiral 
and spin glass phases may belong to different universal- 
ity classes. Furthermore, in contrast to spin glasses, the 
chiral glass phase is likely to exhibit a one-step-like pe- 
culiar replica symmetry breaking [197]. Similar behav- 
ior was recently observed in the chiral glass state of the 
three-dimensional Heisenberg systems [198]. The one- 
step replica symmetry breaking was also reported to take 
place in the random field model [199] and fragile glasses 
[200]. 

In this review we focus on the chiral glass phase in the 
model (4.5) with screening which captures the PME in 
ceramic superconductors. Moreover, the screening effect 
could be substantial in intergranular ordering of these 
materials since the length unit to be compared with the 
penetration depth is the grain size (~ l/xm) rather than 
the short coherence length of the Cooper pair. As the 
screening effect makes the otherwise long-ranged interac- 
tion between vortices short ranged and destabilizes the 
vortex glass (or gauge glass) phase of type-II supercon- 
ductors in a field [201,202]. one wonders if it would even- 
tually wash out a sharp phase transition to the chiral 
glass phase. 



FIG. 27. The temperature and size dependence of the 

Binder ratio of the chirality, gcG, for L = 1. Inset is a mag- 
nified view around the transition temperature Tea — 0.29. 
After Kawamura and Li [27]. 

Because of a rugged energy landscape the multi-loop 
model (4.5) is very hard to equilibrate. This difficulty can 
be partially ovcrcomcd by using, e.g., the Monte Carlo 
exchange method proposed by Hukushima and Nemoto 
[203]. The work of these authors was influenced by the 
simulated tempering method [204], which in turn can be 
understood as a special case of the method of expanded 
ensembles [205]. In the extended ensemble method of 
Rcf. [203] one simulates the sample with a given bond 
realization at Nt distinct temperatures at a time dis- 
tributed in the range \Tmini Tmax] and the whole config- 
urations at two neighboring temperatures arc occasion- 
ally exchanged. Monte Carlo updating consists of the 
two parts: The first part is the standard local Metropo- 
lis updating at each temperature, and the second part 
is an exchange of the whole lattices at two neighboring 
temperatures. Parameters Nt, Tmin and Tmax should be 
chosen in such a way that each replica could wander over 
the whole temperature region and the rate of exchanges 
between two neighboring temperatures would be of order 
of 0.5 and nearly constant [203]. 

We present mainly the results obtained for the induc- 
tance L = 1 by the replica exchange Monte Carlo [203] 
simulations. Fig. 27 displays the size and temperature 
dependence of the Binder ratio of the chirality, gcG ■ The 
data of gcG for I = 3,4,6,8 all cross at almost the same 
temperature T « 0.29, strongly suggesting the occur- 
rence of a finite-temperature chiral-glass transition at 
TcG ~ 0.29 (temperature T is measured in units of J). 
In particular, the data below Tcg show a rather clear 
fan out. The determined value of Tgg is slightly lower 
than the corresponding chiral-glass transition tempera- 
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ture of the pure ±J XY spin glass determined in Ref. 
[196], TcG = 0.32 ± 0.01 (an estimate in [197] yields a 
higher value for Tcg)- Note that the XY spin- glass model 
corresponds to the Z — > limit of the present model. 



The observed suppression of Tcg by the screening effect 
seems reasonable, since the latter makes the long-ranged 
interaction between chiralities short-ranged making the 
chiral-glass transition less favorable. 



TABLE 3. Critical exponents of three-dimensional (3D) Ising spin glass (SG) and of 3D XY chiral glass (CG) with 
and without screening. The standard scaling relations have been used to obtain the full set of exponents from the 
values reported in the original rcfcrencics. 



type 


distribution 


Ref. 


P 


7 


V 


'7 


3D Ising SG 




±J 


[206] 


« 0.55 


« 4.0 


1.7(3) 


-0.35(5) 


3D Ising SG 




±J 


[207] 


« 0.65 


4.1(5) 


1.8(2) 


-0.26(4) 


3D Ising SG 




Gaussian 


[208] 


w 0.64 


w 4.7 


2.0(1.5) 


-0.36(6) 


3D XY CG, L = 





±J 


[196] 


« 0.45 


w 3.6 


1.5(0.3) 


-0.4(2) 


3D XY CG, L = 





±J 


[197] 


« 0.69 


« 2.2 


1.2(0.2) 


0.15(20) 


3D XY CG, L = 


1 


±J 


[27] 


w 0.5 


w 2.9 


1.3(0.2) 


-0.2(2) 




FIG. 28. Finite-size scaling plot of gca with Tea = 0.286 
and ucG = 1.3. After Kawamura and Li [27] 



Fig. 28 shows the scaling function gccil^^^'^'^ \ T — 
Tcg I) given by Eq. (6.11) with fixing Tcg ~ 0.29. 
The best scaling fit gives the correlation length exponent 
ucG = 1.3 ±0.2 [27]. 

The temperature and size dependence of the chiral- 
glass susceptibility, XCG, defined by Eq. (6.8), are shown 
in Fig. 29a. Finite-size scaling analysis based on the 
relation (6.13) is made with fixing Tcg = 0.286 and 
t^CG = 1.3, yielding the chiral critical-point-decay expo- 
nent 77cG = —0.2 ± 0.2. The resulting finite-size-scaling 
plot is displayed in Fig. 29b. Other exponents can be 
estimated via the standard scaling relations (6.14) and 
(6.16) as 7CG ^ 2.9 and /3cg - 0.5 [177]. 



FIG. 29. (a) The temperature and size dependence of the 
chiral-glass susceptibility, Xcg, for L = 1. (b) Finite-size 
scaling plot of xcg with Tcg = 0.286, vcg = 1.3 and 
r]cG = -0.2. After Ref. [177] 

As one can see from Table 3, for L = 1 the exponents 
reported in Ref. [27] are comparable to those obtained in 
Ref. [196] for L = but they are clearly different from 
newer estimates of Kawamura and Li [197]. Therefore, it 
is not clear whether the screening effect is irrelevant or 
not at the 3D chiral-glass transition. Further studies are 
required to clarify this point. 

Within error bars the exponents of chiral glasses ob- 
tained in Ref. [196] and [27] coincide with those for 3D 
Ising spin glasses, suggesting that they belong to the 
same universality class. The results of [197] showed that 
the chiral glass transition may lie in a universality class 
different from the Ising glass. Regardless to this con- 
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trovcrsy, the simulations [27] gave a strong evidence of 
the existence of a chiral glass as a new phase in ceramic 
superconductors . 



n.n - 
-ti.4 
-o.a 

-_.2 
-1.6 



J L 



1 = 6 

1- 1--:' 

J I I I L 



0.0 o.a D.+ o.e o.a 



0.000 



V- 0.001 



-n fin? - 



-0.003 




(T-TJL= 



FIG. 30. (a) The temperature and size dependence of the 
nonlinear susceptibility X2 for Z = 1. An arrow represents the 
location of the chiral-glass transition point, (b) Finite-size 
scaling plot of ;y2 with Tea = 0.286, ucg = 1-3 and 72 = 4.4. 
After Kawainura and Li [177] 

From the experimental point of view the nonlinear sus- 
ceptibility X2 plays a more important role than the chiral 
glass susceptibility XCG (6.10) because experiments can 
probe it directly but not xcG- On general theoretical 
grounds, X2 is expected to show a negative divergence at 
the transition point where the time-reversal symmetry is 
spontaneously broken in a spatially random manner [26] . 
The results shown in Fig. 30a are fully consistent with this 
expectation. The exponent associated with the negative 
divergence is estimated via a finite-size scaling analysis 
(see Fig. 30b) is equal to 72 — 4.4 assuming Tqg = 0.286 
and fcG=l-3. This value of 72 is larger than the chiral- 
glass susceptibility exponent, 7cg — 2.9. At present, it 
is not entirely clear whether this deviation reflects a true 
difference in the asymptotic critical behavior. 

The tendency that the chiral-glass ordering is sup- 
pressed at larger inductances was seen from simulations 
[177]. The phase diagram in the T — L plane is sketched 
in Fig.31. There appears to be a finite critical value of 
the inductance, Lc, above which there is no equilibrium 
chiral-glass transition. Although it is difficult to pre- 
cisely locate Lc due to the extremely slow relaxations at 
low temperatures, it appears to lie around 5 < Lc < 7. 
It was also shown [177] that the paramagnetic tendency 
is more enhanced for larger L, while the chiral-glass or- 
dering itself is suppressed for larger L. 

An attempt to observe the flux glass (or orbital glass) 



was made in Ref. [177] by the Monte Carlo simulations. 
Naively, one expects that the flux defined by Eq. (4.5), 
should behave in the same way as the chirality, since it is 
also a pseudoscalar variable sharing the same symmetry 
as the chirality. Indeed, the flux-glass susceptibility, e.g., 
shows a divergent behavior [177] similar to its chiral glass 
counterpart xcG- However, in contrast to the naive ex- 
pectation, clear crossing of the Binder ratio as observed 
in gcG is not observed in the flux glass case at least in 
the range of lattice sizes I < 10. Rather, the ordering 
tendency seems more enhanced in the sense that the flux 
glass susceptibility gpG tends to increase with increasing 
L exhibiting a feature of the ordered phase even above 
TcG ~ 0.29 [177]. As the flux is an induced quantity 
generated by the finite inductance effect, we believe this 
behavior to be a finite-size effect. Presumably, for in- 
ductance L = 1 the flux hardly reaches its asymptotic 
critical behavior in rather small lattices studied there. 
We tend to believe that the flux or orbital glass can oc- 
cur in the multi-loop model (4.5) of larger system sizes 
and it should accompany with the PME. 
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FIG. 31. A phase diagram in the T-L plane. Renormalized 
inductance L is defined by Eq.(4.3). After Kawamura and Li 
[27] 



C. Differences between chiral glass and other glassy 
phcises in disordered type II superconductors 



The point quenched randomness destabilizes the 
Abrikosov flux line lattice in a pure type-II superconduc- 
tors, yielding novel glassy phases such as Bragg, gauge 
and vortex glass (the so called Bose glass [209] caused by 
columnar or planar defects will not be discussed here). 
The aim of the present chapter is to make a clear dis- 
tinction between these phases and the chiral glass. We 
consider the weak and strong disorder cases separately. 
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TEMPERATURE 

FIG. 32. A phase diagram of a weakly disordered type-II 
superconductor. In three dimensions the vortex glass phase 
should disappear in the presence of the screening. 

Bragg and vortex glass. The starting point to describe the 
Bragg [210-212] and vortex glass [213] is the Ginzburg- 
Landau free energy (5.8) with a weak disorder. The later 
can be introduced via the mean field transition temper- 
ature Tco in the expression of a{T) in Eq. (5.8), i.e, 
Tco ^ TcQ + STcQ (r) , where STcq (r) is a random quantity. 
In the absence of an external magnetic field and disor- 
der the model (5.8) may be mapped to the uniform XY 
model with the specific heat exponent a < [214,215]. 
Then, according to the Harris criterion [216], the weak 
randomness does not change the zero-field critical behav- 
ior of the typc-II superconductors. 

The situation becomes quite different if the external 
field is applied to a sample. In this case the question 
about the H — T phase diagram is still under debate 
but its currently popular view [212,217-219] is shown in 
Fig. 32. As shown first by Larkin [220], disorder spoils 
the translational long-range order in the Abrikosov lat- 
tice since the randomness in the local values of the crit- 
ical temperature leads to a random potential acting on 
the vortices. Although the disordered averaged order pa- 
rameter [< &q){iQ.u) >] (<5 and u denotes a reciprocal 
lattice vector of the Abrikosov lattice and the vortex dis- 
placement, respectively) vanishes in the thermodynamic 
limit, the correlation function 

5(4 r) = [< exp{iQ.[tI(r) - w(0)]} >] (6.17) 

may still obey an algebraic decay [210,217] provided the 
external field is weak enough. Since this algebraic de- 
cay should show up in Bragg peaks the corresponding 



phase is called Bragg glass [211] or elastic vortex glass 
[221]). So the Bragg glass is characterized by transla- 
tional correlations (6.17) decay asymptotically as power 
laws. It reflects some residual order of the Abrikosov 
lattice. After the Bragg glass was first proposed [211] 
(one should be noted that the algebraic decay of S{Q, r) 
of three-dimensional vortex lattices in impure supercon- 
ductors was first predicted by Nattermann [222]) its exis- 
tence was supported by further analytical [223-225] and 
numerical calculations [226,227]. The strong experimen- 
tal evidence for the existence of the Bragg glass phase 
was provided by the neutron diffraction data of Klein et 
al [228]. 

On increasing external field the Bragg glass becomes 
unstable to the vortex glass phase in which translational 
correlations (6.17) decay exponentially. In the vortex 
glass phase, although there is no off-diagonal long-range 
order one may expect a nonzero Edwards- Anderson order 
parameter [213] 

<lv^ = ^J dr-\^{r-)\\ (6.18) 

where gauge-invariant ^{r) = 'l'(r) exp[J^ ?'y4.fiZ] and V 
is the volume of the system. This means that the vor- 
tices in this state arc completely frozen at random po- 
sitions dictated l)y the disorder as spins or chiralities in 
the spin or chiral glass phase. The pinning of vortex lines 
prevents the system from the energy dissipation and the 
linear resistivity should, therefore, vanish in the vortex 
glass phase. Interesting transport properties at and be- 
low the critical temperature were intensively discussed by 
Fisher et al [221] by scaling arguments. Their predictions 
were confirmed by measurements of the I — V character- 
istics [229,230]. The existence of the vortex glass phase 
was also supported by simulations [227] . The evidence of 
the Bragg glass-vortex glass transition was provided by 
experiments [231,232] and simulations [227]. This transi- 
tion appears to be very sensitive to the disorder strength 
[232]. 

Recently, it has been proposed that low- Tc materials 
would be better captured by a phase diagram different 
from what proposed in Fig. 32 [233] . Namely, instead of 
the vortex glass phase the so called multi-domain glass 
would intervene between Bragg glass and disordered liq- 
uid phase. Translational correlations in this new phase 
are argued to exhibit a significant degree of short-range 
order. 

Nonomura and Hu have found the vortex slush phase 
[234] to exist between the vortex glass and disordered 
liquid phases by simulations [235]. Such a new phase 
was not detected by similar simulations of Olsson and 
Teitel [236] but it seems to be supported by the resistiv- 
ity experiments [237,238] and by the analytical argument 
[239]. Further studies are required to to confirm or reject 
the existence of the multi-domain glass and vortex slush 
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phases. 

Gauge glass. So far the disorder was assumed to be weak. 

In the case of strong disorder gauge glass hkc models were 
proposed to describe granular superconductors [240,241]. 
Assuming that each superconducting grain i is described 
by the phase of the order parameter and the Joscphson 
couplings Jij are positive and the same for all of grain 
pairs {Jij = J > 0), the Hamiltonian is given by [202] 

H = -J ^ cos(0i - - A^j - Xo^Uij) + 

<ij> 

Ij^i^xaf. (6.19) 

p 

Here J2<ij> ^^'^ J2p denotes the summation over nearest 
neighbours and over plaquettes of a cubic lattice, respec- 
tively. The influence of the external field and randomness 
is accounted by gauges Aij which are taken to be inde- 
pendent random variables with a distribution between 
and 2n. a^j = jZ' a{f)df, where a denotes the fluctu- 
ations of the vector potential which are limited by the 
bare screening length Aq. 

The gauge glass phase is characterized by a vanishing 
averaged condensate wave function [< >] = 0, but a 
finite Edwards- Anderson order parameter [| < \1' > j'^]. 
The transition to this phase is accompanied with a di- 
verging gauge-glass susceptibility defined in Ref. [188]. 

In the absence of the screening (Aq oo) the gauge 
glass was shown to be ordered at finite temperatures in 
three dimensions, while in two dimensions the ordering 
takes place only at T = [189,190,242-247]. The anal- 
ysis based on the determination of the scaling behavior 
of the domain wall energy £dw{1) oc gave exponent 
9{iD) « [243,242,246], indicating that the 3D gauge 
glass is more marginal than the Ising spin glass which 
have 6{?>D) w 0.2. A detailed discussion of the rela- 
tion between the vortex glass and the gauge glasses was 
givc;n by Blatter et al [248]. The vortex glass occurs in 
the model with weak disorder but the disorder becomes 
effectively strong as the field is increased [249]. If one 
assumes (6.18) as a definition of the vortex glass order 
parameter then the gauge glass would be a vortex glass. 

In ceramic superconductors the length unit associated 
with the intergranular ordering is the moan grain size 
which is of order micron, screening effect is generally non- 
negligible. The influence of screening (finite Ao) on the 
nature of ordering of the gauge glass was considered by 
several groups [201,202,250]. It turns out that screening 
destroys the gauge glass transition in three dimensions. 
This conclusion remains unchanged if one includes the 
effect of anisotropy into the gauge glass model [251,252] 
by assuming an extra contribution to Aij arising from 
the external field. 

Contrary to the Bragg and vortex glass, the gauge 
and chiral glass may occur in systems with strong dis- 



order. Comparing Hamiltonian (4.5) for the chiral glass 
^^ext _ Q-j g^j^i^ (6.19) for the gauge glass one can see two 
main differences. First, the former is supposed to ap- 
pear in the zero external magnetic field, while the latter 
cannot exist without the field which serves as a source 
of frustration. Second, opposed to the gauge glass case 
where the couplings are always " ferromagnetic" , the cou- 
plings in the chiral glass model may be positive or neg- 
ative. This means that the chiral glass may be observed 
only in a system with anisotropic pairing. Since the U{1) 
gauge symmetry is not broken, in this phase, unlike in 
the vortex or gauge glass, the phase of the condensate is 
not ordered, even randomly, on sufficient long length and 
time scales: The thermodynamic ordering occurs only in 
the loop-supercurrents, or in the chiralities. Therefore, 
the chiral glass is not thermodynamically superconduct- 
ing state. The modest screening destroys the gauge glass 
in three dimensions but leaving the chiral glass stable 
against thermal fluctuations. 

D. Experimental search for the chiral glass phase 

First we will discuss some requirements for the ap- 
propriate ceramic samples where one could expect the 
chiral-glass phase. One important parameter character- 
izing the sample is the dimensionless inductance, L, given 
by Eq.(4.3). From Fig. 31 it follows that an equilibrium 
chiral-glass state could be realized in the samples with 
not large L. If one models a loop as a cylinder of radius 
r and height h, its inductance is given by _L = An'^r'^ /h. 
Putting r ^ 1/im, h/r ^ 0.01 and J ~ 20K (these values 
are chosen to mimic the sample used in Ref. [29]), one 
gets L ~ 10~^. Since this value is considerably smaller 
than Lc, an equilibrium chiral-glass phase may well occur 
in such samples. By contrast, if the sample has too large 
a grain size or too strong Josephson coupling, an equi- 
librium chiral-glass phase will not be realized, or at least 
largely suppressed. Another requirement for the sample 
is that the grains must be connected via weak links into 
an infinite cluster, not decomposed into finite clusters. 
This means that the concentration of tt junctions should 
exceed the percolation threshold. The finite-cluster sam- 
ples cannot exhibit a chiral-glass transition, although the 
paramagnetic Meissner effect is still possible [9,253,254]. 

Once appropriate samples could be prepared, the 
chiral-glass transition is detectable in principle via the 
standard magnetic measurements by looking for a neg- 
ative divergence of nonlinear susceptibility X2 [27] as in 
the case of spin glasses. In such magnetic measurements, 
care has to be taken to keep the intensity of applied ac 
and dc fields weak enough, typically much less than IG, 
so that the external flux per loop is sufficiently smaller 
than $0- 

A sharp negatively-divergent anomaly of X2 was re- 
ported in a YB2C4O8 ceramic sample by the ac method 
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by Matsuura et al [29] which might be a signal of the 
chiral-glass transition. This study was extended by 
Deguchi et al [255] to the case when the pressure is ap- 
phcd to a sample to vary the Joscphson couplings. The 
transition temperature Tea was found to be shifted to- 
wards higher temperatures as the pressure increases but 
the character of the transition to the chiral glass state re- 
mains unchanged. This observation is consistent with the 
theoretical prediction [27,177] that Tea is determined en- 
tirely by the strength of the Josephson interaction which 
should grow with the pressure. 

The opposite conclusion from the susceptibility data 
was reported by Ishida et al [32] who showed that the 
nonlinear susceptibility of a (Sro.7Cao.3)o.95Cu02-a; sam- 
ple negatively diverges but it oscillates as a function of 
temperature near the critical point. As argued by these 
authors, such a behavior can be explained by the Bean 
critical state model but not by the chiral glass one. 

Recently, Papadopoulou et al [39] has also observed 
the negative divergence of X2 with exponent 72 ~ 3.4 
for a melt-cast Bi2Sr2CaCu208 sample. This value of 72 
is smaller that the theoretical value 72 ~ 4.4 but close 
to the chiral glass susceptibility exponent jcG ~ 2.9 
[27,177]. The temperature region to observe the diver- 
gence of X2 is, however, too narrow that the result of 
Papadopoulou et al [39] cannot be considered as a good 
evidence for the existence of the chiral glass phase. 

As in the case of spin glasses, measurements of dy- 
namic susceptibilities such as x'{^) x"{^) would also 
give useful information on the possible chiral-glass order- 
ing, particularly when combined with the dynamic scal- 
ing analysis. Near the chiral-glass transition point, the 
imaginary part of the linear susceptibility, x"(w), is ex- 
pected to satisfy the dynamic scaling relation of the form, 

X''(a;,r,if)«.c./^oo/.co.co^.(_|_, ^^J^), (6.20) 

where t =\ {T — Tcq)/Tcg \ and zqg is a dynamical 
chiral-glass exponent. From Ref. [27] (see also Table 3), 
we get the static chiral-glass exponents to be vcg — 1-3, 
/3cG — 0.5 and Pcg + 7CG — 3.4. Using the dynamical 
exponent zcg = 6.3 ± 1.7 [256] we obtain take a value 
around zcC^CG ~ 8.2. 

A dynamic scaling analysis was made by Leylekian, 
Ocio and Hammann for LSCO ceramic samples [257,258]. 
These authors performed both the ac susceptibility and 
the noise measurements, and found an intergranular co- 
operative transition even in zero field at a temperature 
about 10% below the superconducting transition tem- 
perature of the grain. Note that the noise measurements 
enable one to probe truly zero-field phenomena where 
one can be free from the extrinsic pinning effects such as 
the ones envisaged in the so-called critical-state model 
[259]. It was then found that the data of x" satisfied 
the dynamic scaling relation (6.20). Here note that one 



is not allowed to invoke the standard vortex-glass sce- 
nario to explain such intergranular cooperative transition 
in zero field, since in the standard vortex-glass picture 
frustration is possible only under finite external fields. 
By contrast, the experiment seems consistent with the 
chiral-glass picture. 

Meanwhile, when the intragramilar superconducting 
transition and the intergranular transition take place at 
mutually close temperatures as in Ref. [257,258], the 
Josephson coupling, J, which has been assiimcd to be 
temperature independent in the present model, is actu- 
ally strongly temperature dependent in the transition re- 
gion. In such a case, care has to be taken in analyzing 
the experimental data, since the temperature dependence 
of J might modify the apparent exponent value from the 
true asymptotic value to some effective value. In fact, the 
dynamical exponent zv ~ 30 determined by Leylekian 
et al were different from the standard spin-glass value, 
which might be due to the proximity effect of the intra- 
granular superconducting transition [257,258]. If, on the 
other hand, the temperature dependence of J was taken 
into account in the fit, a more realistic value zj/ ~ 10— 15 
was obtained in Ref. [257,258]. As claimed by Leylekian 
et al they data are compatible with the transition to 
the gauge glass state. Using i/qo = 1.3 ± 0.04 and 
zgg — 4.7 ± 0.7 from Ref. [243], we have zcg^gg ~ 6.2. 
The experimental value seems to be more consistent with 
the theoretical estimate zcgvcg ~ 8.2 for the chiral glass 
than with the gauge glass (or vortex glass) counterpart. 
In view of large error bars, such a conclusion remains , 
however, ambiguous. 

The most convincing evidence for the chiral glass phase 
from the dynamical scaling analysis of susceptibility was 
reported by Deguchi et al for YBCO ceramic [260]. Their 
results shown in Fig. 33 give zcg^gg = 8.2 and j3cG — 
0.5 which are well accounted by the chiral glass model. 

Experimental data from the out-of-phase component 
x"{^) of the ac susceptibility and the scaling analysis 
may be employed to obtain the relaxation time which 
diverges on approaching Tcg- In this case, a "freezing" 
temperature Tj, can be defined as a frequency-dependent 
temperatiire Tf{Lo) at which the maximum relaxation 
time is l/w. If there is a phase transition at Tcgi the 
relaxation time t ^ 1/lo behaves as 
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Tfi'^) - Tcg 



Tcg 



t~ 



(6.21) 



as T ^ Tcg- Here Tf{uj) is chosen as the temperature at 
which the peak of x"(a;) is located. With the help of the 
scaling relation (6.21) Papadopoulou et al have obtained 
zcG^CG ~ 11 for the melt-cast Bi2212 sample supporting 
the transition to the chiral glass phase [31]. Moreover, a 
similar analysis on the sintered Bi2212 sample, that does 
not display magnetic aging, gave unphysical results for tq 
and zcG^CG [31]. This result strengthens the evidence 
for the existence of the chiral glass phase. 



31 



' 1 

■ Hi 

:■ 1 IJDI) 



50 



(a) 



60 



10-^ 

I0-' 

](;.' 

10-^ 



H-Ci [GJ 

/i=o.':)i [GJ 

Z i- / 



' 1000 
- 100 



ID" 



10' 



10^' 



FIG. 33. (a) Temperature dependence of the linear suscep- 
tibility for different values of frequency indicated next to the 
curves, (b) The scaling plot for x" which gives zv = and 
P = 0.5. After Deguchi et al [260]. 

Another way to probe the existence of the chiral glass 
phase is based on transport properties [256,30]. As men- 
tioned above, in this phase a sample should not be true 
superconductor, with a small but nonvanishing linear re- 
sistivity PL even at and below TcG- Rough estimates of 
the residual pL were given in [177]. The dynamical scal- 
ing theory for transport properties near the chiral glass 
transition may be developed with the help of its vortex 
glass counterpart [221]. Assume that under the external 
current of density j, there occurs a voltage drop, or elec- 
tric field of intensity E. The voltage drop comes from 
two nearly independent sources: One from the motion of 
integer vortex lines, Ey, and the other from the motion 
of chiral domain walls, E^ [256]. Then we have 



E{j,T) = E,{j,T)+E,{j,T). 



(6.22) 



The first part is essentially a regular part, while the chi- 
ral part should obey the dynamic scaling law associated 
with the chiral-glass transition, i.e., 



(6.23) 



where he spatial dimension has been set equal to D = 3. 
Employing the Fisher-Fisher-Huse argument [221] yields 
the following asymptotic behaviors of the scaling function 



ax, 

a'exp[-6a;~''°G] 



t > 0, 
t < 0, 



as a; — > 0, 

(6.24) 



where a, a' and b are positive constants. /itcG is an un- 
known exponent describing the chirality dynamics in the 
chiral-glass state but it is expected to be positive (in the 
vortex glass case < p, < 1). 

The linear resistivity pL and nonlinear resistivity p^^ 
can be written as a sum of the two nearly independent 
contributions [251], 



PL 



dE 



= PL,v + PL,K , 



3=0 



Pnl 



1 d^E 
6 "df 



PNL,v + PNL,k- 



(6.25) 



j=0 



At the chiral-glass transition point, the vortex terms 
Pl,v{T) and Pnl,v{T) stay finite without prominent 
anomaly, while the chiral parts exhibit a singular behav- 
ior associated with the transition. Then from Eq. (6.23) 
we have 



PlAT) 
PnlAT) f 



t < 0, 

f{zcG-l)-^CG ^ t > 0, 

0, t<Q, 



(6.26) 



Since (zcg ~ ^)'^CG ~ 6.9 is a large positive number 
Pl,k{T) given by Eq. (6.26) vanishes toward Tea sharply 
and the behavior of p^ is dominated by the regular vortex 
term. 

As shown in the next section, Yamao et al have ob- 
served that the nonlinear resistivity and the nonlinear 
magnetic susceptibility of YB2CU4O8 diverge at the tem- 
perature where the magnetic remanencc sets in [30]. 
Meanwhile, the linear resistivity remains finite (see Fig. 
44 below) without any appreciable anomaly. This behav- 
ior is hard to understand if one regards the transition to 
be a transition to the Meissner or the vortex glass phase, 
where pL should vanish. The experimental data obtained 
for the linear resistivity [30] seem to be compatible with 
the chiral glass picture where the total linear resistivity 
remains finite at T = TcG- 

As follows from Eq. (6.26), the nonlinear resistivity 
Pnl shows, however, stronger anomaly than and it 
becomes divergent provided zcg < 5. Using the phase 
representation Kawamura [256] obtained zcg = 6.3 ± 1.7 
for the self-inductance L = 1 and it is hard to justify that 
ZCG < 5. Thus, it is not clear if the nonlinear suscep- 
tibility diverges at the chiral glass transition or not. So 
far as zcg generally depends on the kind of dynamics. 
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it may take a smaller value in the vortex representation 
than in the phase representation. Another possibility is 
that zcG may depend on the self-inductance and other 
values of L could lead to smaller values of zcg- More ac- 
curate estimates of the dynamical exponent are needed 
to make a full comparison with experiments. 

The chiral glass phase is expected to display off- 
equilibrium phenomena similar to those in spin glasses, 
such as aging and memory effects. As mentioned above, 
recently Papadopoulou et al [24] have observed the aging 
effect in ceramic BSCCO sample at very weak fields (see 
Fig. 40 below). It may have some relation to the chiral 
glass order. 

Taken together, the problem of the existence of the 
chiral glass phase in ceramic materials remains open. 
Further theoretical and experimental studies of off- 
equilibrium dynamical and transport properties of the 
chiral glass ordered state would be of much interest. 

VII. DYNAMICAL PHENOMENA: 
EXPERIMENTS AND SIMULATIONS 

In this section we present results of experimental stud- 
ies of the dynamical phenomena reflated to the PME. 
Monte Carlo and Langevin dynamics simulations show 
that all of these phenomena may be captured by simple 
multi-loop model (4.5). 

A. AC susceptibility 

Experiments. It is well known that dc susceptibility pro- 
vides important information on static magnetic proper- 
ties of a investigated material whereas ac-susceptibility 
may shed light on both static and dynamic behaviors 
and provide the clue to the underlying mechanism for 
the PME. 

Fig. 34 shows the ac susceptibility at different frequen- 
cies of a sintered Bi-2212 PME-sample [261] measiucd in 
an ac field Hac — 0.12 Oe at different frequencies. The in- 
phase component x'(w) shows shielding properties that 
directly correspond to the ZFC dc susceptibility. The 
out-of-phase component x"i'^) measures the dissipation 
in the sample and its sharp appearance takes place at 
a temperature very close but somewhat below Tc- The 
onset temperature for dissipation in the sample is al- 
most frequency-dependent within the resolution of the 
measurement [261], only at an appreciably lower temper- 
ature an observable frequency dependence occurs. The 
broad frequency-dependent maximum in x"{^) from Fig. 
34 and the corresponding frequency dependence of x'{^) 
presumably come from spontaneous magnetic moments 
with widely distributed relaxation times. 

In Fig. 35 x' and x" of the above sintered Bi-2212 
sample is shown at a frequency 17 Hz but in different 



superposed dc fields. One remarkable observation here 
is that at temperature close to Tg « 87K the x" curves 
are strongly suppressed even by a dc-field as small as 
0.1 Oe. This behavior contradicts the picture of conven- 
tional flux penetration where an increasing of H^c should 
increase the dissipation and shifts x" to lower tempera- 
tures. The present suppression of the out-phase com- 
ponent may be understood assuming that the low field 
dynamics are governed by the response from the orbital 
magnetic moments. On one hand, the spontaneous mo- 
ments with a relaxation time co, where u) is the angular 
frequency of the ac field, are flipped back and forth by 
the ac-ficld causing the energy loss. On the other hand, 
the dc-field polarizes many of these moments along its di- 
rection preventing them from flipping, and thus reducing 

x". 
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FIG. 34. Ac-susccptibility of a sintered Bi-2212 PME sam- 
ple (Tc — 87 K) at different frequencies: 17 mHz, 170 mHz, 
1.7 Hz, 170 Hz and 1.7 kHz. Prom Magnusson et al [261]. 

The increase of x" at high dc field again is due to 
ordinary intragranular flux penetration [261]. An inter- 
esting difference between sintered and melt-cast samples 
was observed [39,261,262]: the dynamics of the former 
is dominated by spontaneous moments whereas the mo- 
tion of thermally excited vortices plays also an important 
role for the later ones. Thus, the ac susceptibility mea- 
surement points to the existence of spontaneous moments 
responsible for the PME but it does not clarify their ori- 
gin. 

Another way to probe the flux dynamics is to study 
the flux noise spectra S^{f) [263]. The suppression of 
the flux noise in PME samples by the dc field [263,39] 
(see upper panel of Fig. 36) also confirms the existence 
of the spontaneous supercurrents. For large enough dc 



33 



fields the noise is expected to be dominated by vortices 
starting to grow with increasing field. The dc field de- 
pendence of the flux noise resembles that, at the same 
temperature and field ranges, of x". 

In a region where the response of the ac susceptibil- 
ity is linear in field, %" and S^{f) are related by the 
fluctuation-dissipation theorem 



S^if,T,Hdc) = 2kBT 



x"{ f,T,Hdc) 



(7.1) 



where ks is Boltzmann constant. Fig. 36 shows a 
comparison between the magnetic noise and the ac- 
susceptibility via the fluctuation dissipation theorem. It 
is remarkable to note the almost perfect mapping be- 
tween the ac susceptibility and the zero field noise was 
achieved. The validity of the fiuctuation dissipation the- 
orem constitutes that the noise measurements have been 
done in the equilibrium. 
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FIG. 35. Ac-susceptibility of a sintered Bi-2212 PME sam- 
ple a frequency 17 Hz in diflFerent superposed dc-fields. From 
Magnusson at al. [261]. 

Since x"{.f) depends on / weakly, from Eq. (7.1) we 
have 5*0 (/) ~ 1//. The experiments of Magnusson et al 
show that this quantity scales with / as 1//" where a is 
close to 1 indicating that the noise may be a fiicker noise 
and arose from the distribution of activation energies of 
the vortex hopping. Since the 1// noise has been also 
observed in spin glasses [264] one may suggests that the 
PME is relevant to some glassy phase like a chiral glass 
one. 

Simulation. AC susceptibilities may be monitored by ap- 



plying the external magnetic ac field. Then the external 
field H containing the dc and ac parts is given by 



H = Hdc + Hac cos{ujt) 



(7.2) 



In general, the dc field is is necessary to generate even 
harmonics. The real and imaginary parts of n-th order 
harmonics x'ni'^) Xni'^) calculated as 

1 

X'ni'^) = —r— / m{t) cos{nut)d{ut) , 

TT/lac J-7T 



m{t) sin{nojt)d{(jjt) . (7.3) 



For the system described by the multiloop model (4.5) 
the magnetization m{t) is given by Eq. (4.10), the di- 
mensionless ac field hac is related to Hac via (4.12). 
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FIG. 36. (a) Magnetic flux noise fS<s>/T (left-hand scale) at 
0.2 Hz (diamonds), 2 Hz (squares), and 20 Hz (circles) and x" 
at three different frequencies, from left to right 0.17, 1.7 and 
17 Hz (right-hand scale, solid linos) versus T. Filled symbols 
correspond to Hdc = and open symbols to Hdc = 0.3. PME 
sample, (b) Magnetic noise S<s>{H)/S<s>{0) (left-hand scale, 
solid squares) and x" (right-hand scale) at 17 Hz and T = 80K 
vs. Hdc for the PME sample. The maximum at higher fields 
arises from intragranular flux penetration. 5'<i>(7?)/S<i>(0) at 
/ = 17 Hz and T = 62 K versus the superimposed dc field 
for the non-PME sample (solid triangles) is also shown. The 
perfect agreement with the fluctuation dissipation theorem is 
seen. From Magnusson et al [263]. 

The linear ac susceptibility and higher harmonics of 
ceramic s-wave superconductors have been intensively 
studied using the Bean model [265], Kim model [266] 
and a network of resistively shunted Josephson junctions 
[267,268]. Very little theoretical and simulation work has 
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been done in this direction for rf-wave superconductors, 
although such work would be vital to interpret experi- 
mental data [261,263,262]. 

The ac susceptibility of a rf-wave ceramic superconduc- 
tor was computed [145] by Monte Carlo simulations with 
the help of Hamiltonian (4.5) and Eq. (7.3) with n = 1. 
While Monte Carlo simulations involve no real dynam- 
ics, one can still expect that they give useful information 
on the long-time behavior of the system. In fact, the 
characteristic time for the sintered samples, which are 
believed to be captured by our model, is of order 10~^^s 
[267]. This time has the same order of magnitude as a 
single Monte Carlo step. In general, the period of oscil- 
lations is much longer than the characteristic time [145]. 
For such a slowly changing ac field the system can be 
regarded as being in quasi-equilibrium and the Monte 
Carlo updating may be applied. A priori, the validity of 
this approximation is not clear but it may be justified 
[269] by comparing the Monte Carlo results with those 
obtained by other approaches to the dynamics such as 
considered in Ref. [267]. 

A qualitative agreement between simulations [145] and 
experiments [261,262] was achieved, e.g., for the depen- 
dence of x' S'lid x" on temperature and the ac field. 
Namely, as Hac is increased, %" grows and its peak shifts 
towards lower temperatures, while %' becomes less and 
less diamagnetic. Simulations [145] showed that at high 
frequencies there is very little difference between the x of 
the s- and d-wave models. At low frequencies, clear differ- 
ences have been seen: The magnitude of x" of the s-wave 
model becomes an order of magnitude smaller than that 
of the d-wave model, indicating that the d-wave system 
exhibits much stronger dissipation. A close similarity 
observed between the x" of d-wave superconductors and 
that of the spin glass [145] is consistent with the chiral 
glass picture as well as with experimental findings [261]. 

It would be highly desirable to check the validity of the 
fluctuation dissipation theorem (7.1) and verify whether 
the 1 //-noise occurs in a system showing the PME by 
simulations. This may be done in the framework of the 
multiloop model (4.5). 



B. Compensation effect 



Experiment. Heinzel et al. [23] have shown that the PME 
may be analysed by the compensation technique based on 
the measurement of the second harmonics of the mag- 
netic ac susceptibility. The essence of the so called com- 
pensation effect is that the second harmonics signal van- 
ishes at some compensation field Hcom even a nonzero 
dc field is applied. This is a bit surprising because even 
harmonics should be generated in the presence of the dc 
field due to M{H) = -M{-H) symmetry. 
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FIG. 37. Second harmonics of the magnetic ac susceptibil- 
ity {Hac=l Oe) obtained for YBa2Cu307 after field cooling 
in a dc field Hdc=30 Oe to T = 77K and stepwise reducing 
Hdc- (a) Modulus \x2\ of the complex harmonics, (b) Real 
(X2) and imaginary (^2) parts of the second harmonics. After 
Heinzel et al. [23]. 



The compensation effect may be detected in the fol- 
lowing way. The sample is cooled in the external dc field 
down to a low temperature and then the field is switched 
off. At the fixed low T the second harmonics are moni- 
tored by applying the dc and ac fields to the sample. 

The results obtained for PME sample YBa2Cu307 are 
shown in Fig. 37 where Hcom ~ 15 Oe. The observed 
intersection of x'2 ^^'^ X2 Hdc — Hcom indicates a 
change of their sign. The key observation is that the 
compensation effect appears only in the samples which 
show the PME but not in those which do not. It may 
be explained in the following way. Due to the presence 
of non-zero spontaneous orbital moments and the frus- 
tration, the remanent magnetization or, equivalently, the 
internal field is generated in the cooling process. If the 
direction of the external dc field is identical to that dur- 
ing the FC procedure, the induced shielding currents will 
reduce the remanence. Consequently, the absolute; value 
of the second harmonics |x2| decreases until the signal of 
the second harmonics is minimized at a field H^c = Hcom ■ 
Thus the study of the compensation effect allows us to es- 
tablish not only the existence of the spontaneous orbital 
moments but also to determine internal fields generated 
during the cooling process. The compensation effect is 
a collective phenomenon and it may be captured by the 
multi-loop model of and tt-J unction network. 
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FIG. 38. The second harmonics of the d-wave model (4.5) 
obtained after field cooling in a dc field hdc ~ f from T = 0.7 
to T = 0.1. At the lowest T = 0.1 (< Tea ^ 0.17 for Z = 4) 
the dc field used in cooling is switched off and the second 
harmonics are generated by applying the combined dc and 
ac fields. The dc field is stepwise reduced from hdc = 1 to 
hdc ~ 0. The inductance is chosen to bo equal to L = 4 
and system size I = 8. The arrows indicate the sense of the 
changes in the dc field. The results axe qualitatively the same 
as those presented in Fig. 37. After Ref. [269]. 

Simulation. We now present the Monte Carlo simulation 
results [269] explaining the compensation effect [23] in 
ceramic superconductors. Again we employ model (4.5) 
and Eq. (7.3) (n = 2) to compute the second harmonics. 
The calculations follow exactly the experimental proce- 
dure of Hcinzcl et al [23]. In the FC regime, first the 
system is cooled in some dc field hdc 7^ from a high T 
down to a low T which is below the paramagnet-chiral 
glass transition temperature [269] . When the lowest tem- 
perature is reached the dc field used in cooling is switched 
off and we apply the combined field given by Eq. (7.2). 
We monitor the second harmonics reducing the dc field 
from a high dc field to zero stepwise. Fig. 38 shows the 
typical behavior of the second harmonics for the d-wave 
ceramic superconductors. \x2\ reaches minimum at the 
compensation field hcom = 0.7 ± 0.05. At this point, 
similar to the experimental findings presented in Fig. 37 
[23], the intersection of x'2 ^^'^ X2 is observed. This fact 
indicates that at Hcom the system is really in the com- 
pensated state. Furthermore, in accord with the experi- 
ments, at the compensation point the real and imaginary 
parts should change their sign [23] . Our results show that 
x'2 changes its sign roughly at hdc = hcom- A similar be- 
havior is also displayed by X2 but it is harder to observe 
due to a smaller amplitude of X2- 



Performing simulations in the same way as for the d- 
wave case, one cannot observe the compensation effect 
for the s-wave model [269] because the cooling of a un- 
frustrated " ferromagnetic" system could not produce any 
remanent dc field which would compensate the external 
one. This result is again in accord with the experimental 
data [23]. 

Thus, the compensation effect may be explained, at 
least qualitatively, by using the multiloop model of the 
ceramic superconductors with anisotropic pairing sym- 
metry. The existence of the compensation phenomenon 
only in those samples showing the PME indirectly sup- 
port the chiral glass phase in these materials. 



C. Aging phenomenon 



ExperimeMt. The aging phenomenon observed first in 
spin glasses [270] have been studied both theoretically 
[271] and experimentally [264] in detail. In this phe- 
nomenon the physical quantities depend not only on the 
observation time but also on the waiting time (or age), 
tw, i.e. how long one waits at constant field or temper- 
ature before measurements. The origin of such memory 
phenomena relates to the rugged energy landscape which 
appears due to disorder and frustration [144]. 

Overall, one can explain the aging effect qualitatively 
in the following way. After the waiting time twi and tw2, 
e.g., the system gets trapped , in principle, in different lo- 
cal minima as shown in Fig. 39. Since the corresponding 
energy barriers are different the physics should depend on 
in what local minimum the system was before one starts 
to do measurements. It should be stressed that both 
the waiting and measurement times are typically much 
shorter than the relaxation time which may be astronom- 
ically large for a frustrated system [144]. The aging is, 
therefore, a nonequilibrium phenomenon. 

Recently, Papadopoulou et al [24,272] has observed the 
aging effect in the PME sample Bi2Sr2CaCu208 monitor- 
ing the ZFC magnetization. The relaxation of the ZFC 
magnetization has been measured by cooling the sample 
in zero field to the measuring temperature, allowing the 
sample to stay at that temperature for a certain time tw 
and then applying the probing field and recording the 
change of the magnetization with observation time at 
constant temperature. Papadopoulou et al have made 
two key observations. First, the aging effect is not ob- 
served at high fields and at high temperatures. Second, 
this effect occurs at some intermediate temperature in- 
terval (see Fig. 40) and it disappears again at low tem- 
peratures. Thus the aging phenomenon exists only for 
weak enough fields and in an intermediate temperature 
region. 
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FIG. 39. Schematic energy landscape of a frustrated sys- 
tem in the phase space. After waiting times twi and tw2 the 
system gets trapped in different local minima. 

The first observation is trivial because at high fields 

or high temperatures the roughness of the energy land- 
scape does not play any crucial role and the system looses 
its memory. At low temperatures the role of the energy 
landscape becomes important, the second result of Ref. 
[24] is, therefore, not trivial from the point of view of the 
standard spin glass theory. Papadopoulou et al suggested 
that at low temperatures the external field is screened 
from the bulk of the sample and it cannot probe the col- 
lective behavior of the Josephson junction network. As 
one can see below, the correctness of this idea may be 
confirmed using the multi-loop model (4.5). 

The results from aging measurements are in immedi- 
ate contrast to models explaining the PME as a flux com- 
pression phenomenon [19,21] which are only applicable to 
the positive FC magnetization observed in some conven- 
tional superconductors. The aging phenomenon in the 
melt-cast samples [24], with similarities to the aging in 
spin glasses, implies the interaction between spontaneous 
magnetic moments caiising them to behave collectively. 
Inside each of the large grains of a sample there can ex- 
ist many thousands of domains, and the intragrain re- 
gions can be, therefore, modeled as Josephson networks 
containing a random distribution of ordinary Josephson 
junctions and tt junctions. 

Simulation. The study of the aging phenomenon in spin 
glasses by computer simulations lasts for may years [273] . 
Our aim is not to present those results but to focus on 
the explanation [274] of specific features of the aging in 
granular superconductors observed by Papadopoulou et 
al [24]. 

To understand the experiments [24] we use the model 
(4.5) for a d-wave superconductor. The dependence of 
the magnetization m{t) given by Eq. (4.10) on the wait- 
ing time was computed by the Monte Carlo simulations 
[274]. In order to mimic the aging effect in the ZFC 
regime we quench the system from a high temperature 



to the working temperature. The system is there evolved 
in zero field during a waiting time, tw Then the exter- 
nal field h is turned on and the subsequent growth of the 
magnetization M{t,tw) is monitored. 
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FIG. 40. M^fc/Hdc of the melt-cast Bi2Sr2CaCu2 08 sam- 
ple is plotted vs observation time for Hdc = 0.02G at tem- 
peratures (a) 65 K, (b) 82 K, and (c) 86 K. The waiting time 
used were 300 and 3000 s; in (b) the result for twait = 30000 
s is also included. Prom Papadopoulou et al [24]. 

One can demonstrate [274] that there are three screen- 
ing regimes for the aging phenomenon. In the strong 
screening limit when L > L\ = 9±0.5, the aging is 
suppressed at any temperature. In the weak screening 
regime, Z < = 3.5 ± 0.5, it is observable even at 
low temperatures. The intermediate screening regime, 
L* < L < Lj, is found to be the most interesting: the 
aging is present only in an intermediate temperature in- 
terval and it does not appear at low temperatures. 

Fig. 41 shows the results for Z=7 belonging to the 
intermediate screening regime. In agreement with the 
experiments [24] , the aging effect appears only for the in- 
termediate temperature interval. At low T's [T <T* = 
0.02 ± 0.01) the effect is suppressed due to the screen- 
ing of the magnetic field from the bulk. The results have 
been obtained for the observation times comparable with 
the waiting ones but we believe that they should be valid 
for longer observation time scales [274]. 

It should be stressed that the experimental finding of 
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Papadopoulou et al [24] cannot be explained by the stan- 
dard XY model where the screening effect is not taken 
into account. The mechanism of aging may be under- 
stood from studies the spatial distribution of flux inside 
a sample [274]. In fact, in the strong screening limit the 
aging effect does not occur at any temperature because 
the external field is screened entirely. For a fixed screen- 
ing strength the magnetic field is expelled more and more 
from the bulk as T is lowered. Therefore, for the interme- 
diate values of screening one could not observe the aging 
at low T's. 

The fact that the increase of the self-inductance would 
wash out the aging effect may be explained by explor- 
ing the effect of screening on the energy landscape [275] . 
Studies of local minima at T = showed that the energy 
landscape gets smoother and smoother as the screening 
is enhanced and the glassy effects would become less pro- 
nounced. 
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FIG. 41. The temperature and time dependence of M for 
t„ = 1000 and = 25000. Z = 24, L = 7 (intermediate 

screening regime) and h = 0.05. The aging disappears at low 
temperatures. The results are averaged over 60 - 120 samples. 
Prom Li et al [274]. 



So, the non-trivial aging phenomenon in ceramic super- 
conductors may be understood in the framework of the 
model (4.5) where the anisotropy of the pairing function 
plays a key role. The observation of aging and the PME 
in the same material does yield support, therefore, for 
the existence of d-wave superconductivity. Furthermore 
the PME dynamics owns the characteristics features of 
the non-equilibrium dynamics of spin glasses. The re- 
sults on aging are thus indicative of the formation of a 
low temperature glassy phase in a system showing a pro- 
nounced PME, a phase that resembles the chiral glass 



phase [26,27]. 



D. Anomalous microwave absorption 

Experiment. MWA measurements are a sensitive tool for 
detection of the transport properties of high- Tc super- 
conductors below the critical temperature, when super- 
conducting paths are risen and dc measurements cease 
to give information. Soon after the discovery of high- 
Tc superconductivity it has been established [276-280], 
that in ceramic samples there exists a large non-resonant 
MWA in zero and low applied magnetic field that can 
be associated with the transition to the superconducting 
phase. In fields below Hd the MWA of high-Tc materials 
is due to dissipative flux motion in the link network or, 
equivalently, due to the junction's resistance. 

Braunish et al. [3,4,281] found a striking correlation be- 
tween the field dependence of the MWA and the PME. 
Their typical results are shown in Fig. 42. Clearly, in the 
absence of the PME the MWA has a conventional mini- 
mum at H = and in low fields P{H) ^ H^. The change 
over from the quadratic to linear behavior appears at H*^ 
associated with the critical field, at which the magnetic 
flux breaks into the weak- link network [4]. Contrary to 
non-PME samples the MWA of PME samples show a 
maximum at i? = and it reaches minimum at a field 
near which the FC signal changes from paramagnetic to 
diamagnetic behavior [282]. 
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FIG. 42. (a) The field dependence of the MWA for a pow- 
dered HTSC sample which does not show the PME. The 
changeover from parabolic to linear occurs around 0,6 Oe. 
(b) The same as in (a) but for the PME sample. The minima 
are located at ±0.6 Oe, near the field where the FC signal 
changes from diamagnetic to paramagnetic, (c) The MWA 
signal of a Bi-2223 sample showing a weak PME. From Ref. 
[4]. 

Simulation. The unusual behavior of the MWA may be 
explained in the framework of a single- loop model [9]. 
If the loop contains a tt junction, then the spontaneous 
internal field is generated. In the absence of the exter- 
nal field, the microwave field is coupled to this field and 
dissipates energy. As a static field is turned on, the os- 
cillatory field loses gradually its access to the sponta- 
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neous moments and the absorption power decreases. If 
the dc field is strong enough the tt junction loop would 
behave like a junction one, leading to the increase of 
the power. The absorption minimum should be located 
at H « $0/45* [9], i.e. between the fields making the 
centers between the first maximum and minimum of the 
free energy shown in Fig. 19. 

The general scheme to study the MWA in an Joseph- 
son junction networks may be found, for example, in Ref. 
[283]. Here we follow a simpler approach [284]. Namely, 
in order to reproduce the anomalous behavior of the ab- 
sorption power we use the resistively shunted junction 
(RSJ) model [285] combined with the linear response the- 
ory [286]. In this formalism the linear response to the 
external electromagnetic field is proportional to a volt- 
age - voltage correlation function. To get insight on the 
MWA we calculate quantity P which is proportional to 
the frequency integral of this correlation function. Then, 
we have [284] 



< Vf > 



(7.4) 



where < > is a mean value of the square of the volt- 
age induced by the thermal noise on each junction and 
R is the normal resistance of the links. 
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FIG. 43. The field dependence of P for s- (left panel) and 
d-wave (right panel) ceramic superconductors. We choose 
T = 0.2 and L = 0.1, 1 and 10. The results are averaged 
over 20 samples. After Li [275]. 

To calculate Vi we use the RSJ model [12] for the cur- 
rent flowing between two grains. The derivation of dy- 
namical equations for this model is given in Appendix B. 
The full set of equations for the gauge invariant phases 9 
(see Eq. (B.2) is as follows 



2eR dt 



2e 

— J^(n) sin 61^(11) 



-—A- [A+0,{n) - A+0.(n)] - Un,t) . (7.5) 

Here we use following notation: the site of each grain is 
at position n = (n^;, n^, n^) (i.e. i = n); the lattice direc- 
tions are n = x,y, z; the link variables are between sites 
n and n + /i (i.e. link ij = link n, /i); and the plaquettes 
p are defined by the site n and the normal direction /x (i.e 
plaquette p = plaquette n, /x, for example the plaquette 
n, z is centered at position n -|- (x -I- y)/2). The forward 
difference operator A+0i,(n) = 9^{n + ^) — 0^{n) and the 
backward operator A~^i,(n) = 6v{n) — 6^{n— ^) (see Ap- 
pendix B for more details). Critical currents Jp(n) are 
the same as in Hamiltonian (4.5). The left part of Eq. 
(7.5) describes the normal current, while the first and 
second terms in the right part correspond to the super- 
current and screened current, respectively. The Langevin 
noise current (n{n,t) has Gaussian correlations 



(C^(n,i)C^,(n',i')) = 



.'Sn,n'S{t-t') . (7.6) 



The local voltage Vi is then given by 

dt 



(7.7) 



Eq. (7.5) describes the overdamped dynamics. One 
can show that the inertia (capacitive) terms do not 
change results qualitatively and they are neglected. 

The system of differential equations (7.5) is integrated 
numerically by a second order Runge-Kutta-Helfand- 
Greenside algorithm for stochastic differential equations 
[287]. The time step is chosen to depend on L and is 
equal to At = O.lrj and At = O.lrj x Z for Z > 1 and 
L < 1, respectively. 

The field dependence of the MWA in two-dimensional 
s- and d-wave disordered superconductors was studied 
by Dominguez et al [284]. Similar results have been ob- 
tained for three dimensions [275] . Fig. 43 shows the field 
dependence of the MWA for T=0.2 and for various values 
of L. In the case of s-wave superconductors we have the 
standard minimum at H = for any value of inductance 
and T. As (expected, P ~ at weak fields. 

For the d-wave samples P has the unconventional peak 
a,t H = 0. Contrary to the one-loop model [9] such peak 
is seen not only for L > 1 but also for L < 1. 

In our model (7.5) the temperature dependence of the 
critical current is neglected. However, one can show that 
the dimensionless temperature T chosen in Fig. 43 corre- 
sponds to the relevant to experiments real temperature, 
Tji. In fact, the critical current depends not only on 
temperature but also on conditions under which sam- 
ples were prepared. The typical value of the critical cur- 
rent density for ceramic superconductors is ^ 10^ A/ 
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( see, for example, Ref. [280,288-290]). Since the typical 
size of grains is about l/xm we have the critical current 
7c ~ W~^A. Using Tr = JT/ks = hIcT/2ekB one ob- 
tains Tr/T ^ lOQK . Clearly, the dimensionless T chosen 
in Fig. 43 correctly describes the experimental values of 
temperature [3]. 

Comparing Fig. 43 with Fig. 42 one can sec that 
the multi-loop model (4.5) correctly captures the exper- 
iments of Braunish et al [3,4]. Furthermore, the height 
of the peak shown in Fig. 43 is very small: (P{H = 
0) — Pmin)/Pmin is of Order of 10""^. This is also in qual- 
itative agreement with experimental findings [3] that the 
peak should be low. 



E. AC resistivity 

Experiment. In recent experiments Yamao et al. [30] have 
measured the ac linear resistivity, po, and the nonlinear 
resistivity, p2, of ceramic superconductor YBa2Cu408. 
po and p2 are defined as the first and third coefiicient of 
the expansion of the voltage V{t) in terms of the external 
current Iext{t): 



The experimental value of the power law exponent is 
a« 1.1. 



V = poIext + P2lext + ■■ 



(7.8) 



When the sample is driven by an ac current Iext{t) = 
Iq sin(cji), one can relate po and p2 to the first harmonics 
and third harmonics V^^ in the following way 

= V{t)sm{n^t)d{wt) , n = 1,3 . (7.9) 

Fig. 44 shows the temperature dependence of V^/I 
for YBa2Cu408 sample around Tc2 below which the re- 
manent magnetization becomes nonzero [30]. The key 
observation from this figure is that po does not vanish 
even at and below the intergrain ordering temperature 
Tc2. On the other hand, p2 has a peak near this tem- 
perature, which was found to be negative [30] (see Fig. 
45). 

In the chiral glass phase the U{\) gauge symmetry is 
not broken and the phase of the condensate remains dis- 
ordered. [27] The chiral glass phase, therefore, should not 
be superconducting but exhibit an Ohmic behavior with 
a finite small resistance. Based on these theoretical pre- 
dictions Yamao et al. [30] speculated that their results 
give further support to the existence of the chiral glass 
phase, in addition to previous results from magnetic sus- 
ceptibility measurements [29]. 

Another interesting result of Yamao et al. [30] is the 
power law dependence of \Vl^{Tp) / Iq)^\ (or of P2) at its 
maximum position Tp on Iq: 
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FIG. 44. The temperature dependence oiV^/ 1 around Tc2. 
Frequency / = 23 Hz and values of / are shown in the figure. 
The inset shows the data measured at 7 = 0.05 mA. The solid 
line indicates the Umiting po — T curve. After Yamao et al 
[30]. 



Simulation. In order to reproduce the experimental re- 
sults of Yamao et al. [30] we again use the multi-loop 
model (4.5) combined with the Langevin dynamics equa- 
tions of the RSJ model which are given by Eq. (7.5). We 
should add to the right part of these equations the ex- 
ternal oscillatory electric field Eext{t) cx lext [291]. The 
details of simulations are given in [291]. 



Fig. 46 shows the temperature dependence of the lin- 
ear resistivity po = V^/Iq and the non- linear resistivity 
P2 = -4:V^^^{T)/I§ for different values of Iq. The linear 
resistivity is nonzero at the chiral glass transition temper- 
ature Tea as well as at the point Tp where the nonlinear 
resistivity has maximum. The results for po presented in 
Fig. 46 were obtained for the system size I = 8 but one 
can show [291] that they are valid in the thermodynamic 
limit. Thus, the fact that po does not vanish at Tp is 
in full agreement with the experiments [30]. The other 
observation from simulations [291] is that Tp coincides 
with a temperature where the FC magnetization starts 
to departure from the ZFC one. This is also consistent 
with the the experimental findings of Yamao et al. [30]. 
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FIG. 45. The temperature dependence of V^^^/I^. Fre- 
quency / = 23 Hz and values of / are shown in the figure. 
The peak is located at Tp w 27.5 K. After Ref. [30]. 

As said above Yamao et al. have made an interest- 
ing conjecture that the peak of the nonlinear resistivity 
should correspond to the phase transition to the chiral 
glass phase. Tea and Tp are, however, well separated 
(see Fig. 46) suggesting that the peak in the nonlinear 
resistivity may not be indicative for the existence of the 
chiral glass state. One of possible interpretations of sim- 
ulation results is that Tp just separates the normal state 
phase from a "chiral paramagnet" where there exist local 
chiral magnetic moments. These moments can be polar- 
ized under an external magnetic field, an therefore one 
can observe the PME under a low external field below 
Tp. At a lower temperature, collective phenomena due 
to the interactions among the chiral moments will start 
to be important, leading to the transition to the chiral 
glass state. This last transition should show up in the 
nonlinear chiral glass susceptibility which diverges at Teg 
[145]. 

At present, it is not entirely clear whether Tp corre- 
sponds to the chiral glass transition. Its departure from 
TcG may be just a shortcoming of the RSJ dynamics. 

From lower panel of Fig. 46 one can see that the height 
of the peak in max\V^^^/lQ \ grows as the current ampli- 
tude decreases supporting the experimental result pre- 
sented in Fig. 45. A detailed study [292] shows that this 
dependence can be described by the power law (7.10) 
which is valid not only for both d-wave and s-wave ce- 
ramic superconductors. The exponent a was found to 
be not universal but depend on the self-inductance and 
current regimes [292]. 

In the weak current regime a is independent of the 
self-inductance and a = 0.5 ± 0.1 for both of s- and d- 
wave materials. In the s-wave case, since Tp corresponds 
to a true continuous metal - superconductor phase tran- 



sition, one can use the scaling argument [221] to obtain 
the dynamical exponent z via a [292]. In three dimen- 
sions 2; = 5 — 2a Ri 4 which is higher than z ^ 3 fov 
the corresponding uniform model [202]. For the rf-wave 
system, Tp is the temperature where there is an onset of 
positive magnetization, i.e. the paramagnetic Meissner 
effect starts to be observed, but it docs not seem to cor- 
respond to a phase transition [291]. The scaling analysis 
is not, therefore, applied. 
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FIG. 46. The temperature dependence of V^/Io and 
Vsui/Io of a d-wave superconductor for I = 8, L = 1 and 
u) = 0.001. The open triangles, squares and hexagons cor- 
respond to lo = 0.1,0.05,0.02. The arrows correspond to 
Tp = 0.8 and Tec = 0.286, respectively. The results are 
averaged over 15 - 40 disorder realizations. Prom Li and 
Dominguez [291]. 

In the strong current regime [292] the exponent a de- 
pends on the screening strength. For the d-wave super- 
conductors and 1 < i < 5 we have awl which is close 
to the experimental value [30] . In order to make the com- 
parison with experiments to be meaningful we show that 
a was measured in the strong current regime. In fact, the 
real current is / = %^/o, J ~ 10^ K and Iq ~ 10~^ we, 
therefore, have / ~ 10"^ mA which is much weaker than 
the current / ~ 10 mA used in experiments of Yamao et 
al. The another issue is that the interval of inductance 
where a « 1 is realistic for ceramics [280] because typical 
values of L are bigger than 3 [280] . An accurate compar- 
ison between theory and experiments requires, however, 
the knowledge of L which is not known for the compound 
of YBa2Cu408 studied in Ref. [30]. 



F. Enhancement of critical current 

Experiment. The question about enhancement of the 
critical current of superconductors is important from 
both the theoretical and application points of view. Sev- 
eral experiments [293,294] in ultrathin (5 - 10 nm thick) 
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YBa2Cu307-rE films showed that an electric field E may 
affect Tc as well as the critical current Ic- The later may 
increase or decrease depending on the polarity of E and 
this is usually attributed to the changes in the critical 
current due to variations in the charge density or due to 
a redistribution of carriers which appear at the surface 
layer with depths of the order of the electrostatic screen- 
ing length dE (in high-Tg superconductors, cIe ~ 5^)- 
For any experimentally measurable field- induced effects, 
dE should be larger than the superconducting coherence 
length and this happens to the the case in ceramics which 
have a reduced carrier density [293-295]. 

Recent experiments [25,296-302] revealed, however, 
that the enhancement of Jc is also possible in bulk (1.5 
mm thick) ceramic high-Tc superconductors when E is 
applied trough an insulating layer (see Fig. (47). This 
unusual behavior cannot be explained by a surface effect 
[293,294]. Moreover, for high enough electric fields, the 
critical current always increases regardless of the polarity 
of the field. 
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FIG. 47. I-V characteristics for (BiPb)22Sr2Ca2Cu30^ 
sample for S = (1) and E = 120 mV/m (2). J^, = 320 
A/cm^ and T = 77 K. After applying the electric field the 
IV curve is shifted towards higher currents, i.e. the apparent 
critical current increases. After Smirnov et al [300]. 



Using the SNS contact model one can show that an 
electric field can induce a change in the critical currents of 
the Josephson junctions in granular samples [303] . How- 
ever, in this model there is either an enhancement or a 
deleption of depending on the sign of E. Authors of 
Ref. [304] have proposed that an electric field applied to a 
granular superconductor can produce a magnetoelectric- 



likc effect, which could be indirectly related to the behav- 
ior of the critical current observed in [296-299,305-307], 
but no comparison with the experiments was given. 
Simulation. Following Domingucz et al [308] one can ex- 
plain the experimental findings on enhancement of the 
critical current of bulk granular superconductors using 
the rf-wave m\ilti-loop model. In fact, we could employ 
the full Langevin equations (7.5) but the unusual electric 
field effects may be captured even neglecting the screen- 
ing effect {L = 0) and thermal fluctuations (77;^(n, t) =0). 
Then the current J/j(n) between to grains n and n +jl 
is a sum of the Josephson supercurrent plus a dissipative 
Ohmic current [308] 



J^(n) = J°_^sin6'^(n) 



^0 de^jn) 
2ttR dt ' 



(7.11) 



where ^/^(n) = + jl) — (j)(r\) — An{n,t) is the gauge in- 
variant phase difference, and A^(n, t) — /""'^'^A.dl. 
Together with the conditions of current conservation, 
^^[J^(n) — J^(n — /I) = 0, equations (7.11) defines dy- 
namics of the system. 

Wc assume that an electric field E is applied in the z 
direction and the network is driven by an external cur- 
rent density I^xt along the y direction. Then, in the ab- 
sence of the external magnetic field, the vector potential 
A^(n, t) = —S^^ziJEt — 5fi^yay{t), where the electric field 
frequency loe = 2TrEd/^o with d being the intergrain 
distance or junction thickness [308]. As follows from Eq. 
(7.11), the external current density determines the evo- 
lution of ay{t) as 



lext = -^^Iny 0y{n) + 



$0 day 
2'kR dt 



(7.12) 



The average voltage per junction induced by the driving 
current is then obtained asV = (<I>o/27r) < day/dt >. 

It should be noted that one docs not take into ac- 
count the effects of intergrain and intragrain capacitances 
in model (7.12) and the screening of the electric field 
is, therefore, neglected [308]. The electric field scale is 
Eq ^ RIo{T = 0)/d = 7rA(0)/2e«i, where A(0) is the 
superconducting energy gap: for the YBaCuO ceramics 
one has A(0) « 20 meV, d ^ 10 - 20l, which gives 
Ef) ~ 30 MV/m, i.e. in the same range as the fields used 
in [296-299,305-307]. Dominguez et al [308] studied two 
models of disorder, (i) Granular s-wave superconduc- 
tor: is assumed to be a random variable uniformly 
distributed in the interval [/o(l — Ac),/o(l + Ac)] with 
< I^^ >= lo and Ac < 1. (ii) Granular d-wave supercon- 
ductor involving and tt junctions has the same bimodal 
distribution for as discussed in previous sections. 

Fig. 48 shows typical simulation results [308] for 
current- voltage characteristics before and after applying 
an electric field E for a s— (left panel) and rf-wave (right 
panel) three-dimensional ceramic superconductor. In the 
former case after the field is switched on the whole IV 
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curve shifts to lower values of the current indicating the 
decrease of the apparent critical current and the voltage 
change AV = V{E) — V{0) is positive for a given cur- 
rent I > Ic- This contradicts, however, the experimental 
finding shown in Fig. 47, where an increase in Ic was 
seen. 
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FIG. 48. Current-voltage characteristics before and after 
applying an electric field E. (a) For a granular s-wave super- 
conductor with Ac = 0.6. (b) For a d-wave superconductor 
with a bimodal distribution of positive and negative Joseph- 
son couplings coming from and tt junctions. Voltages are 
normalized by I Rio and currents by l^Io and the lattice size 
I = 16. After Dominguez et al. [308]. 

In the d-wave case, after application of the electric field 
the IV curve moves towards higher currents, i.e. the ap- 
parent critical current increases. It is precisely what was 
observed on experiments (see Fig. 47). Furthermore, if 
one changes E —E the IV curves overlap showing that 
the effect is independent of the polarity of the electric 
field. The similar effect was also seen in the experiments 
[298,299]. 

There is other interesting aspect of experimental find- 
ings that if an applied current is near the critical current 
then AV/Vb increases as a function of E and, after reach- 
ing a maximum, it falls and then becomes negative for 
large enough fields [297]. This behavior can be also ex- 
plained in the framework of the d-wave model [12]. 

Two unexpected results from the experiments on the 
electric field effects are: (i) a pronounced electric field 
effect in a bulk ceramic sample and (ii) an increase of the 
apparent critical current as a function of E independent 
of the field polarity. All these features are qualitatively 
captured by the simple d-wave model in which the frus- 
tration due to TT junction and the ac Josephson effect 
induced by the electric field are important. It would be 
interesting to explore the effect of the electric field on the 
paramagnetic signal, on the chiral glass phase and glassy 
dynamics in ceramic superconductors [298,299]. 



VIII. SUMMARY 

We have reviewed the experimental facts, simulations 
and theories of the PME observed in conventional and 
ceramic high- Tc superconductors. There exist two main 
mechanisms for this very interesting phenomenon: the 
flux compression and the tt junction. The later probably 
occiirs due to anisotropic pairing of electrons but other 
possibilities like the scattering on magnetic impurities in 
Josephson junctions are not excluded. 

From the technical point of view, the paramagnetic 
signal due to the fiux trapping inside a sample of con- 
fined geometry may be obtained either by considering the 
Bean critical model or by solving the nonlinear Ginzburg- 
Landau equations. In the case of ceramic superconduc- 
tors with the rf-wave symmetry of the order parameter, 
the spontaneous magnetic moment leading to the PME 
occurs in the loop of odd number of tt junctions. The 
screening plays a crucial role in both the single loop and 
multi-loop models. It shoidd be noted that the later can 
capture not only the PME but also many related dynam- 
ical phenomena in granular materials. 

The flux compression scenario should work for meso- 
scopic conventional materials but it may be also applied 
to granular superconductors where the grain surfaces can 
play an important role. Therefore, the existence of the 
PME itself could not serve as an unambiguous indica- 
tor for the rf-wave symmetry of the order parameter of 
cuprates. It would be very important to prepare ce- 
ramic samples with negligible surface effects to check if 
the paramagnetic signal can occur or not. 

In addition to the phase-sensitive and phase-insensitive 
experiments, the experiments on the aging, anomalous 
MWA, 1// flux noise, compensation effect and enhance- 
ment of the critical current which may be described by 
simple XY model with screening (4.5), tend to support 
the rf-wave pairing symmetry of the cuprate superconduc- 
tors. From both experimental and theoretical points of 
view it would be very exciting to find out whether these 
effects appear in conventional superconductors showing 
the PME. 

The chiral glass phase which may be observed by neg- 
ative divergence of the nonlinear susceptibility is an in- 
teresting example of the time-reversal symmetry break- 
ing. In contrast to the gauge and vortex glass, it appears 
only in a zero external magnetic field. Furthermore, the 
moderate screening could not destabilize the chiral glass 
phase at finite temperatures in three dimensions. Since 
the experiments [32,29 31,255,260] gave conflicting re- 
sults on the existence of this phase in high- Tc supercon- 
ductors, much more work in this direction is left for future 
studies. 
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APPENDIX A: CRITICAL SELF-INDUCTANCE 
OF THE MULTI-LOOP MODEL 

In order to demonstrate that the critical inductance 
of the multi-loop model is smaller than that of the 
single-loop model, we analyse the ground state of an 
isolated frustrated plaquette, an elementary unit of our 
lattice model as shown in Fig. 49. The plaquette con- 
sists of four bonds (junctions) among which one is a tt- 
junction with negative coupling — J < 0, and other three 
are 0-junctions with positive coupling constant eJ > 
(0 < e < oo). The multi-loop model considered here 
corresponds to e = 1, while e = oo describes a single tt 
junction analysed by Sigrist and Rice [9]. The external 
field is set equal to 0. 
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FIG. 49. An isolated frustrated loop consisting of four junc- 
tions, one of which is a tt junction of magnitude -J (thick 
line) and the other three are 0-junction of magnitude e J with 
< e < oo. The case e = 1 corresponds to an elementary 
frustrated plaquette of the multi-loop model, while e — oo 
corresponds to the frustrated single-loop model analysed in 
the previous section. 

Introducing gauge-invariant phase difference on each 
bond Oij = (pi — — Aij ("temporal gauge" [12,145]), 
the dimensionless energy of the plaquette shown in Fig. 
49 is 



H/J = € cos 9i — cos 02 — cos — cos 64 + 

^ 92 + 93 + 04)^. (A.l) 



Here, for simplicity, we drop one sub index of the 9ij 
and denote the phases as 9i at tt junction and 02 — 04 at 



jimctions. Minimizing F with respect to 9i — 0^, we 
get four coupled equations. After a little algebra, one can 
see that the ground-state configuration should satisfy the 
relations 



92 = 9s = 94 = 9, 9i = arcsin(— esin^). 



(A.2) 



with 9 being determined from the following equation 

esin6' + i(36i-Farcsin(-esin6i)) = 0. (A.3) 
L 

One can show that if e < 3, this equation always has dou- 
bly degenerate solutions = ±9* which corresponds to 
supercurrents with opposite directions. A spontaneous 
moment, therefore, arises for an arbitrary value of induc- 
tance. If e is larger than three, there appears a finite 
critical value of the inductance I'c(e) = 1 — 3/e [145]: for 
a inductance L > Lc{(), there again exist doubly degen- 
erate solutions 9 = ±0* , while for L < Lc{e), only trivial 
solution exists at ^ = corresponding to the state with- 
out spontaneous supercurrent. In the limit e — > 00, one 
recovers the critical value Lc = 1 obtained by Sigrist and 
Rice [9]. Thus, we have shown that, contrary to the one- 
loop model, the critical inductance of the interacting loop 
models with weak links of nearly the same magnitude is 
equal to Lc = 0. 



APPENDIX B: RSJ MODEL 



In this appendix we derive the equations of the RSJ 
model which have been used to study dynamical phenom- 
ena such as the MWA, the ac resistivity and the enhance- 
ment of the critical current in chapter VII. As shown in 
the inset of Fig. 50 the RSJ model of a single junction 
contains the superconducting and resistive parts which 
are denoted by a cross and a rectangular, respectively. 
Then the total current I through the junction is [7] 



I ^ Ir + Is 



^^d0_ 



L sin. 9, 



(B.l) 



where R and /<; are the resistivity and the critical current 
of the junction, respectively. The gauge invariant phase 
6 between points 1 and 2 is 

9,2 = ^^1.(v0-^a), (B.2) 

where V(j> is the gradient of the phase of macroscopic 
wave function, A is the vector potential, and the line in- 
tegration is taken between two sides of the junction with 
sign convention in accordance with the coordinate direc- 
tions shown in Fig. 50. 



44 







1 z 






i,j,k: iJ-t ^ 


J ij-l.,L- 

/ 


1 


1 ^ / 


as / 












X 







^ilk ■■ < i < I ~ 1;0 < j < I - 1;0 < k < I; 
^Hk ■■ < i < I - l;0 < j < l;0 < k < I - I; 

:0<z</;0<j<^-l;0<fc<;-l. (B.6) 



The corresponding restrictions on {i,j,k) for phases are 



'rj,fc:0<i<Z-l;0<i<Z;0<fc<Z; 
'fj,fc:0<i</;0<j<Z-l;0<fc<Z; 
?fj,fc :0<i<Z;0<j</;0<fc</-l. (B.7) 



Assuming the self-inductances of all loops to be iden- 
tical and equal to L one has 



ij,k - for (i,j,k) e (B.6) 

otherwise. 



FIG. 50. Schematic description of the three-dimensional 
RSJ model of the Josephson junction network. The gauge in- 
variant phases are marked on some representative junctions. 
Tiio 6'^ values grow in the positive directions the the corre- 
sponding a^xis jj, = x,y,z. The positive direction of currents 
circulating in loops is chosen to be counterclockwise. The 
inset shows a single junction containing the resistive and su- 
perconducting parts. 

We now consider the three-dimensional network on a 
cubic lattice shown in Fig. 50. Representing the current 

flowing from point (i,j,k) in direction fi to the nearest 
point as I^j ^, from Eq. (B.l) we have [285] 

where ^C^- ^ is the gauge invariant phase between points 
{i,j,k) and its nearest neighbor in direction ^. Ifj^k is, 
e.g., the current from k) to {i + 1, j, k) (see Fig. 50) 
and I^- f.^ is the local critical current. 

The crucial point here is that I-^- ^ consists of contribu- 
tions from four loops sharing the corresponding junction 
and we, therefore, have 
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(B.4) 



Labeling the flux ^fj^k threading through loop k) 
in /xi^-plane in the same manner as lij^k ^® have 
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(B.5) 



For the cubic array containing I x I x I junctions fluxes 
(B.5) are defined for following values of k): 



Here <I>^^* is a flux threading through a loop in /xi^-plane 
which is equal to = HS if the /xt^-plane is perpen- 
dicular to H and = otherwise (see also Eq. (4.6)). 
Substituting Eq. (B.4), (B.5) and (B.8) intoEq. (B.3) we 
obtain the full set of equations for gauge invariant phases 
which is often used to study dynamics of the Joseph- 
son junction network at zero temperature. For nonzero 
temperatures one has to add the Langevin noise term to 
mimic the effect of the heat bath. 

We now show that in the case when the external mag- 
netic field H = the Langevin equations for phases may 
be written in a compact form. Using Eq. (4.7) and 
$0 = hc/2e one has 



i,j,kc i,j,k' 



(B.9) 



where J^- is the Josephson couplings. Prom Eqs.(B.3), 
(B.4), (B.'5), (B.8) and (B.9) we obtain the following 
equation for d^j^- 



2eR dt 



'^^[i^lj,k + 



+hj,k 



„;„ nx 
ijk "i,j,k 

nx — ny \ 

' ^i,j+l,k ^i,j,k) 
' ^i,j,k ~ ^i,j-l,k) 



2eL 

+ (^fj,fe-l + ^tj,k ~ ^i+l ,j,/c-l ^ ^'tj,k-l) 
~{(^lj,k + ^i,j,k+l ~ ^i+l,j,k ~ ^i,j,k)-] 



(B.IO) 



Equations for and 6^ j j, can be written in a similar 

form. For the sake of convenience we represent (i,j, k) by 
vector n and redefine 6'^ . ,. ^^{w) and J^- ,. J^(n). 
We introduce also the forward operator and back- 
ward operator A~ which act on ^^(n) in the following 
way 



A+0,(n) =e,(n + ^)-e,(n), 
/^-e,{n)=e,{n)-e,{n- ii). 



(B.ll) 
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In the (i, j, A;)-representation the last equation takes the 
form 

for component ^ — x. Similar equations for /i = y and z 
can be easily written down. With the help of Eq. (B.12) 
one can show that Eq. (B.IO) is reduced to Eq. (7.5) in 
the main text. 
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